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Abstract 

The discovery of a 125 GeV Higgs boson and rising lower bounds on the masses of superpart¬ 
ners have lead to concerns that supersymmetric models are now fine-tuned. Large stop masses, 
required for a 125 GeV Higgs, feed into the electroweak symmetry breaking conditions through 
renormalization group equations forcing one to hne-tune these parameters to obtain the correct 
electroweak vacuum expectation value. Nonetheless, this fine-tuning depends crucially on our as¬ 
sumptions about the supersymmetry breaking scale. At the same time, U{1) extensions provide 
the most compelling solution to the y problem, which is also a naturalness issue, and allow the 
tree-level Higgs mass to be raised substantially above Mz- These very well-motivated supersym¬ 
metric models predict a new Z' boson which could be discovered at the LHG, and the naturalness 
of the model requires that the Z' boson mass should not be too far above the TeV scale. Moreover, 
this fine-tuning appears at the tree level, making it less dependent on assumptions about the su¬ 
persymmetry breaking mechanism. Here we study this fine-tuning for several 17(1) supersymmetric 
extensions of the Standard Model and compare it to the situation in the MSSM where the most 
direct tree-level hne-tuning can be probed through chargino mass limits. We show that future LHG 
Z' searches are extremely important for challenging the most natural scenarios in these models. 
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I. INTRODUCTION 


The discovery of an approximately 125 GeV Higgs mia at the Large Hadron Collider 
(LHC) has interesting implications for physics beyond the Standard Model (SM) and sn- 
persymmetry (SUSY). On the one hand, it provides a light Higgs boson, as expected from 
snpersymmetry, and can be fitted in the minimal snpersymmetric standard model (MSSM). 
On the other hand, the Higgs mass is slightly heavier than the constrained version of the 
MSSM (cMSSM) can accommodate natnrally [31 H]. 

In the MSSM the Higgs mass canses a natnralness problem because at tree level it has 
an upper bound of the mass of the Z boson, Mz- The dominant higher-order corrections 
to the Higgs mass come from stops, and to obtain a 125 GeV Higgs they need to be rather 
heavy. Heavy stops will provide a large contribution to the low-energy value of the soft 
breaking mass for the up-type Higgs scalar, through the evolution of the renormalization 
group equations (RGBs) from the grand unihcation (GUT) scale to the electroweak (EW) 
scale. This affects the SUSY prediction of the electroweak vacuum expectation value (VEV), 
n, or Mz- This naturalness problem motivates both further examination of nonminimal 
SUSY models that can raise the Higgs mass without the need for heavy stops and alternative 
possibilities for how the soft breaking parameters get generated, which might set them at 
lower energies, reducing the influence the stops have on 

In addition to that naturalness issue, often referred to as the little hierarchy problem, 
the MSSM also suffers from the /i problem. This is also a naturalness problem since there 
should be a natural explanation of how the fi superpotential parameter can be set to the 
same scale as the soft breaking masses. 

f/(l) extensions of the MSSM provide a very elegant solution to this fi problem [5HT2] 
and also raise the Higgs mass with new F and D terms. Nonetheless, as was recently 
demonstrated in the context of the exceptional supersymmetric standard model (EgSSM) 
[T3Hr5] . such models can still suffer from naturalness problems with the mass of the new 
Z' associated with the break down of the new U{1) appearing in the electroweak symmetry 
breaking (EWSB) conditions at tree level [16]. Despite this the constrained version of the 
EgSSM (cEgSSM) [T71 [T8| was still found to be significantly less tuned than the cMSSM. 
Tree-level fine-tuning from the Z' mass was also considered previously |19j . 

However, this comparison of fine-tuning depends crucially upon the assumptions of these 
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gravity mediated SUSY breaking motivated constrained models and, in particular, the uni¬ 
versality constraints being applied at the GUT scale. As mentioned above, given the hndings 
at the LHC, it is worth considering other possibilities, which may allow the soft masses to 
be set at lower energies. As the scale at which the parameters fulhll some breaking inspired 
constraints is lowered the stop masses contribute less to the hne-tuning. 

At the same time in U{1) extensions, lowering the UV boundary scale for the RGE 
evolution also allows even larger F-term contributions to the Higgs mass, so long as one 
only requires A, the coupling between the Singlet Higgs, S and the up- and down-type Higgs 
bosons, Hu and Hd, to remain perturbative up to the UV scale and not all of the way up to 
the GUT scale. 

However, the tuning from the Z' mass limit does not disappear as the UV boundary 
condition is lowered. This tuning appears in the EWSB conditions at tree level and is quite 
difficult to avoid without introducing a pure gauge singlet [20] . 

In this paper we investigate how big this tuning is if we bring this scale all the way 
down to 20 TeV, effectively minimizing the contribution from the stops. We find that the 
Z' limit is enough to already require moderate fine-tuning in the EgSSM. We also show 
this is comparable to the situation in the MSSM dehned at the same scale if charginos 
could be ruled out below 700 GeV. We then show how this tuning from the Z' mass looks 
for different U{1) extensions. Ending that the current severity depends upon the charges 
but that Z' limits are important in constraining the most natural scenarios of these models. 
Therefore, the Z' constraint is amongst the most important in terms of tuning and attacking 
natural supersymmetry experimentally and the next run of the LHG will be crucial in this 
respect. 

Finally, we make a case study, for a few benchmarks, of the impact of raising the high- 
scale boundary condition, Mx, at which the SUSY breaking parameters must be fixed by 
some SUSY breaking mechanism. We show that which model has less fine-tuning depends 
on Mx- We also see rather complicated behavior in the tuning for the EgSSM points due to 
the combination of different sources of tuning. 

As mentioned earlier, the-fine tuning of the cEgSSM was recently studied [16] and there 
it was revealed that the associated Z' boson leads to a new source of fine-tuning since its 
mass appears in the EWSB conditions. 

However in this study we will examine this source of fine-tuning in more detail by con- 
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sidering low-energy constructions where the usual hne-tuning problem from the Higgs is 
minimized. We will also consider alternate charges for the extra U{1) symmetry to relax the 
focus on the EgSSM and demonstrate that this is quite a generic result. 

To quantify the hne-tuning we will employ the traditional Barbieri-Giudice measure m 
[22]. This has been used extensively within the literature e.g. Refs. [iniEiaHi. 

A number of alternative measures have also been applied in the literature [501 - 165] with 
varying motivations. A very different approach is to work within a Bayesian analysis. There 
the concept of naturalness is automatically incorporated since in models where one must 
hne-tune parameters to ht measured values of the observables, the region with high likelihood 
will occupy a tiny prior volume |H 1661170] suppressing the posterior. Indeed in the MSSM 
and the next-to-MSSM (NMSSM) if one transforms GUT scale parameters to the VEVs, the 
inverse of the Jacobian for this transformation looks quite like the derivatives that appear 
in the traditional hne-tuning measure [661 1671 170] . If one thinks more generally, then a 
model without hne-tuning is one where the parameterization is such that all the parameters 
are observables [011 ED]. This provides a quite general dehnition of hne-tuning as 1/|J| 
where |J| is the determinant of the Jacobian for the coordinate transformation between 
the parameters and the observables. Interestingly this means the tuning is the ratio of the 
inhnitesimal observable space volume element to the inhnitesimal parameter space element 
and coincides with the measure proposed in Ref. [63] when the interval of variation is taken 
to zero. 

While this approach has many merits here we will employ the traditional measure of 
hne-tuning because it is both simple to apply and easy to compare with previous results due 
to it’s widespread use. Fortunately the derivatives which appear in these tunings are also 
similar to the Bayesian motivated measure so there should not be too large a discrepancy 
between the two approaches. 

The structure of this paper is as follows. In Sec. [^we review the models we consider. In 


Sec. Ill we specify the EWSB conditions of the models, with particular focus on how the Z' 


mass inhuences the prediction of Mz- Then in Sec. IV we introduce our hne-tuning measure 
and our approach to evaluating it to obtain the individual sensitivities. The results are then 
given in Sec. |V| 
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II. U{1) EXTENSIONS AND THE EgSSM 


In this paper we consider U{1) extensions of the MSSM where the gauge group at low 
energies is 

SU{3)c X SU{2)w X U{1)y x (1) 

U{iy is the new gauge group beyond that of the SM and MSSM. The minimal superheld 
content of 17(1) extensions which solve the /i problem should be ordinary left-handed quark 
Qi and lepton Li {i = 1, 2, 3) superhelds along with right-handed superhelds hf, dy, ef {i = 
1, 2, 3) for the up-type (s)quarks, down-type (s)quarks and charged (s)leptons respectively 
and three Higgs superhelds, up-type down-type Hd and a singlet under the SM gauge 

group S. 

Here we will refer to 17(1) extensions of the MSSM, which solve the /i problem, as the 
USSM [HHI2]- The couplings for the 17(1)' gauge group should allow the following renormal- 
izable superpotential terms required in the USSM, 


WussM = ■ Qi + vrAQ, ■ Hd + vmLi ■ Hd + XSHd ■ H^, 


D 


( 2 ) 


with h j e {1, 2, 3}. For the 517(2) dot product we follow the convention A-B = ea/sA'^iS^ = 
- A^BL 

The 17(1)' charges should allow for cancellations of gauge anomalies. The most elegant 
way to do this is to use an extra 17(1) gauge symmetry that can be obtained from the break 
down of the Eq gauge symmetry which is anomaly free and have all matter helds that £11 
the three generations of 27-plet representations of Eq survive down to low energies. Such 
models are often referred to in the literature as Eq inspired, and we will adopt this here. 

The breaking of Eq into 50(10) gives rise to Eq —)■ 50(10) x 17(1)^, and the subse¬ 
quent breaking of 50(10) into 517(5) gives 50(10) —)■ 517(5) x 17(1)^ (this is reviewed in 
e.g. Ref. [ZI])- The extra 17(1) gauge symmetry at low energies should then be a linear 
combination of these in the Eq inspired case. 


17(1)'= [/(!);, cos 0 + 0(1)^ sin 0. (3) 

In Table the charges for several popular Eq inspired 17(1) extensions are shown. 

17(1) and Eq inspired extensions of the MSSM have been studied very widely in the 
literature [211 IT2IIM] (or, for reviews, see Refs. [711195]). There has also been a lot of work 
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TABLE I. The C/(l)y, U{1)^ and U{1)n charges of the chiral superfields in the Eq model. 

The specific case of U{l)]\f, corresponding to the EgSSM, is obtained for 6 = arctanVTs. 

recently including investigations of the neutralino sector [96]-E9]; the relic density of dark 
matter |100j : GUT scale family symmetries which can explain the hierarchy of masses in 
the fermion sector and their associated mixings IMl; neutrino physics HEl; explanations of 
the matter-antimatter asymmetry of the Universe though EW baryogenesis or leptogenesis 
[93t IHU 1103] : decays of the Z' boson |in4flin7] : dipole moments |in8j : anomaly mediated 
SUSY breaking with D-term contributions [109] and the (extended) Higgs sectors |lin[llll] . 

Here we will focus most on the special case where the gauge symmetry is U{1)n, under 
which the right-handed neutrino does not participate in gauge interactions. This is the 
case in the EgSSM [TSHIS], and closely related variants pm irmiTTH] . Since the right-handed 
neutrino has no gauge symmetry protecting it’s mass from becoming extremely heavy such 
models may explain the tiny observed masses of neutrinos via the see-saw mechanism and 
the baryon asymmetry in the Universe via leptogenesis [93l 111711118j . Recently it has also 
been studied in the context of electroweak baryogenesis |119] . 

The gauge coupling running in the EgSSM at the two-loop level leads to unihcation more 
precisely than in the MSSM |12Uj or, in slightly modihed scenarios, two-step unihcation 
can take place |112l I121j . If the exotic particles are light in these models this can open up 
nonstandard decays of the SM-like Higgs boson [201 1122[ I123j . 

The correct relic density could be obtained entirely through an almost decoupled “inert” 
neutralino sector m- However, this is no longer phenomenologically viable due to limits 
from direct detection of dark matter and due to a signihcant suppression of the 

decay of the lightest Higgs boson into SM states, due to a new channel into inert singlinos 
opening up. 
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There are still several remaining options. One may specialize to scenarios known as the 
EZSSM |115] where the inert singlinos that cause these problems are entirely decoupled and 
the relic abundance is fitted with a binolike candidate with a novel mechanism involving 
back-scattering into a heavier inert Higgsino. Another well motivated scenario admits two 
possible dark matter candidates |116] , where one will be an inert singlino and the other will 
have a similar composition to MSSM neutralinos. The simplest phenomenologically viable 
solution in that case is to make the singlinos extremely light hot dark matter candidates, 
in which case the lightest ordinary neutralino accounts for almost all of the observed relic 
abundance. 

The impact of gauge kinetic mixing in the case where both of the extra U{1) symmetries 
appearing from the breakdown of Eq are present at low energy was studied in Ref. |128] . 
The EgSSM was also included in studies looking at how Erst- or second-generation sfermion 
masses can be used to constrain the GUT scale parameters |129] and the renormalization of 
VEVs |13U11131j . The particle spectrum and collider signatures of the cEgSSM have been 
studied in a series of papers, [m OKI Oroni [132]. The threshold corrections to the DR gauge 
and Yukawa couplings in the EgSSM have also been calculated and their numerical impact 
in the constrained version examined |133] . 

With three generations of matter 27-plet representations of Eq surviving to low energies, 
the low-energy matter content in each generation, after integrating out the heavy right- 
handed neutrinos, includes. 


(4, *•, dt, 4, + (A, A) + + m + {Hf), 


( 4 ) 


where the Si, and Hf have the quantum numbers of a SM singlet, and up-, down-type 
Higgs helds, respectively, and the Di and Di are SU{3)c triplets that reside in the same 
SU{5) multiplets as these Higgs-like states. 

If one wishes to maintain gauge coupling unification this set of states should be augmented 
by two extra SU (2) doublet states E[' and h' belonging to other 27' and 27^ multiplets that 
must be incomplete at low energies. 

The full superpotential for Eq inspired models coming from 27 ® 27 0 27 decomposition 
of the fundamental Eq representation will then be 


We, = IUo + fGi + W2, 


( 5 ) 
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where 


= XijtSiHf ■ HI + . L» 

+ygt«f//" ■ Qi + yfJlQt ■ Hf + • Hf, (6) 

H'l = g?iAQi ■ Qk + (7) 

n-'j = gS,N!D,dl + + g,%Qi ■ LjD,. (8) 

Nonetheless, while this model is very elegant so far, the snperpotential of Eq. (§ contains 
dangerons terms which can induce proton decay and lead to large flavor changing neutral 
currents (FCNCs). There are a number of approaches to suppress these terms, involving the 
use of different discrete symmetries. Here for the purposes of renormalization group running 
we will simply include the following unsuppressed snperpotential terms, which follows the 
approach taken in work on the cEgSSM HZIIIS], 

W ^ yrLs • Hdcl + ybQs ■ Hddl + ytHu ■ Q 3 U 3 

+XiSHf ■ + KiSDiti + iiH' ■ W, (9) 

where we denote by = Hd, = Hu and S 3 = S the third-generation Higgs and 
SM singlet helds that are assumed to acquire nonzero VEVs. In addition to the terms 
coming from the 27 ® 27 ® 27 interactions given in Eq. ([^, this snperpotential also contains 
a bilinear term fi'H’ ■ if', arising from 27' ® 27^, which is invariant with respect to the 
low-energy SM gauge group and the additional H(l)' symmetry and also anomaly free. 
This term is responsible for setting the masses of the components of the superhelds i7', H', 
included to ensure gauge coupling unihcation, but it is not involved in the process of EWSB. 
Consequently, the impact on the fine-tuning of the value of /r' is much smaller than that 
coming from other sectors, and so can be safely neglected in our study. In all of the scans 
we present below the value of /i' is hxed to /r' = 5 TeV. 

III. ELECTROWEAK SYMMETRY BREAKING 

The Higgs scalar potential for the Eq models considered can be written as [13] 

V = Vf + Vd + V,^,, + AV , (10) 

where 

Vp = X^\S\WHd\^ + \Hu\^) + X^\Hd ■ Hu 
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( 11 ) 


Vn = ^ {\Hu\‘^ - + ^(gilZ/^r + Q 2 \Hu? + Qs\S\Y: (12) 

o Z Z 

Vsoft = ml\S\^ + + m^jH^ + [XA^SH, ■ + H.c.]. (13) 

In these expressions g 2 , g' = \J?)/bgi, and g'^ are the SU{2), (non-GUT normalized) U{1)y 
and U{iy gauge couplings, respectively, and g"^ = g 2 + g'^. The charges Qi, Q 2 and Qs are 
effective U{1)' charges for Hd-, Hu and S, respectively, and A = A 3 . In the case of the U{1)^ 
model, Vf may also contain an elementary p term, as occurs in the MSSM. The term AV 
contains the Coleman-Weinberg contributions to the effective potential. For the purposes of 
this study, we include in AV only the one-loop contributions from the top quark and stop 
squarks. 



Explicit expressions for the running DR top mass rrit and stop masses ^ given below. 
Demanding that the Higgs fields Hi, H 2 and the singlet S acquire real VEVs of the form 



m ^ ^ 


. (s) ^ 


(15) 


at the physical minimum leads to the minimization conditions 

^ o A^ , , 9 , XAx 9^/9 9x ^ dAV , , 

/l = TUh^Vi + —iv 2 + S^)vi - - -^{vj - vf)vi + Dh^Vi + = 0, (16a) 

/2 = + y (g^ + S^)V2 - ^SVi + |-(t^2 - ^'i)^'2 + DhuV2 + = 0, (16b) 

/s = mis + y (^2 + + Dss + = 0. (16c) 

The quantities Dh^ and Ds appearing above are U{iy D-term contributions that are 
absent in the MSSM and NMSSM and are given by 

Dtj, = y {Qivl + Q 2 V 2 + Qip- (17) 

We also include these U{iy D-term contributions in the diagonalized stop masses. 
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where Xt = At — ^ire soft breaking scalar masses and At is 

a soft trilinear coupling. By definition we take mi^ to correspond to the lighter of the two 
states. 


As was noted in Ref. [16] , the hrst two of the conditions in Eq. (16) may be rewritten in 
the form 






+ 


m 


H, 


m 




tan^ (3 Dh^ - Dh„ tan^ (3 


2 

sin 2/3 = 


+ 


tan^ /3 — 1 

y/2XAxs 


tan^ (3 


™t/. + + Dh^ + Dh. ’ 


(19) 


( 20 ) 


with — vf + V2 and tan/3 — V2/V1 and where we have for convenience 

absorbed the effects of the loop corrections into the soft masses, 

1 dAV 


rriTT = m. 


rrijj = m 


H 2 


+ 


+ 


vi dvi 

1 dAV 
V 2 dv2 


Written in the form of Eq. (19), we see the potential new source of hne-tuning alluded to 
above, in the form of the third term on the right-hand side. For large values of the VEV s, the 
D-term contributions can be quite a bit larger than M|. In particular, recent experimental 
limits |134] require that the Z' mass be large, with for example bounds of Mz' ^ 2.51 TeV 
in U{1)^ models and Mz' > 2.62 TeV in U{1)^ models. To satisfy these limits typically 
requires large values of the singlet VEV s. For example, s > 6 TeV is required in the EgSSM 
with U{iy = U{1)n, so that iD^y, \Dh 2 \ ^ Vf| for Eq models with Qs 7 ^ 0 . As a result 


the remaining terms on the right-hand side of Eq. (19) must be tuned to cancel this very 
large contribution to Mz- Moreover, because this is a large tree-level hne-tuning, it may 
negate the improvement in naturalness that is associated with having a reduced need for 
heavy superpartners. In U{1) extended models for which Qs 7 ^ 0, the importance of the Z' 


mass to the hne-tuning in these models can be made even clearer by writing Eq. (19) in the 
form given in Ref. 


c{9, tan /3) 


Ml 


2„2 




+ 


mrr — rriTT 
Hd ^22 


tan^ (3 — 1 


tan^ (3 ^ ^ Ml, 

- + (i(6',tan/3)——. 


( 21 ) 


where 


c{9, tan (3) = 1 — 


9? 


(tan^ (3 — 1) g 


{Qi - Q 2 tan^ (3) [Qi cos^ (3 + Q 2 sin^ (3) , ( 22 ) 
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d{6, tan /3) 


Qi - Q2tan^ /3 
Qs (tan^ /5 - 1) ’ 


(23) 


Written like so, it is evident that the hne-tuning contribution coming from the new D terms 
depends both on the U{iy charges and the Z' mass, and that the tuning can be expected 
to increase with Mz'- As shall be shown below, the exact size of this tuning then depends 
strongly on the choice of U{iy charges, via the coefficient d. 

The extra U{iy gauge symmetry may mix with the U{1)y gauge symmetry associated 
with hypercharge through gauge kinetic mixing. 


pkin _ 

/I -T* 


c\ fiiy 5 


(24) 


where and are held strengths associated with the U{1)y and respectively. 

The gauge kinetic mixing can have a signihcant impact on the phenomenology |135fll37] 
and may reduce the Z' mass limit. 

However, if the extra U{1) gauge symmetry appears from the breakdown of Fg, then sin y 
should be zero at the GUT scale. Nonetheless, even if this term is zero at the outset, it will 
still be radiatively generated if the trace of the U{1) charges, Q'v nonzero. In the 

cases studied here, the trace of the charges over states in the complete 27-plets vanishes, 
but to be consistent with single-step gauge coupling unihcation, we also included H' and 
H' which lead to a nonzero value for YhiQjQ'i- The value induced by this at the EW scale 
though is rather small as can be seeiQin Fig. 3 of Ref. |135] . and this was also checked with 
two-loop RGBs in the EgSSM [T31IIH]- For this reason and due to the huge expansion in the 
number of terms in the two-loop RGBs when one allows for gauge kinetic mixing, we will 
neglect this in our analysis here and throughout this paper. 

In general, though, it is possible for gauge kinetic mixing to be much larger, which can 
be the case if one considers an additional 5-1-5 pair of SU{5) multiplets |135] or which has 
been looked at in the U (1) b-l |136l 1137] . In such a case, this will impact the results in two 
ways, hrstly by altering the Z' limit from experiment and secondly by altering the charges 


which appear in the EWSB condition, which can be seen from examining Eqs. (21)~(23). 


^The specific incomplete multiplets we consider here correspond to the third of the four possible embeddings 
referred to in Ref. [135] . 
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IV. THE FINE-TUNING MEASURE 


As stated above, to quantify the resulting fine-tuning we apply the traditional Barbieri- 
Giudice measure [211 [22] • A specihc model is characterized by a set of n model parameters 
{pi} and is dehned at some input scale Mx- For a given parameter p in this set, one computes 
an associated sensitivity, 

_ a In M| _ 

^ d\np M| dp 

The coefficient Ap measures the fractional variation in M| resulting from a given variation 
in the parameter p. The overall hne-tuning is then taken to be A = max {Ap.}. 

The sensitivity Ap may be calculated directly from the expression for M| in terms of the 
Pi for a particular model, which leads to a so-called master formula for calculating the hne- 
tuning. A master formula for the EgSSM, obtained from the tree-level scalar potential, was 
presented in Ref. [I6] . In order to derive the expression presented there, the fact that s S> n 
was made use of to neglect certain O(n^) terms in the EWSB conditions, greatly simplifying 
the hnal result. For the purposes of exploring a wider class of Eq inspired models, we have 
derived the master formula without neglecting any 0{v‘^) terms. The more complete tree- 
level master formula is somewhat complicated. This is because, unlike in the MSSM, even at 
tree level it is not possible to solve explicitly for the VEVs vi, V 2 in terms of the Lagrangian 
parameters. It may be written in the form 






dq 

dp 


( 26 ) 


where the sum is over all low-energy running parameters appearing in the tree-level EWSB 
conditions, i.e., q G {X, gi, g 2 , g'l}. Expressions for the quantity C and 

the Aq appearing above are given in Appendix]^ It should be noted that the effects of U (1) 
mixing are neglected in deriving Eq. ( [2^ . 

However, it is well known in the MSSM that radiative corrections can signihcantly change 
(indeed, reduce) the hne-tuning |138] . It is, therefore, important when studying the hne- 
tuning to include loop corrections to the ehective potential in the hne-tuning measure. To 


do so it is most convenient to work in terms of the EWSB conditions Eq. (16), rather than 


Eq. (19). The general procedure that we use is as follows (this method has also previously 
been applied in the NMSSM; see, for example. Ref. |139j h For a model in which m helds 
develop real VEVs (e.g. m = 2 in the MSSM, m = 3 in the NMSSM and in the Eg models 
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considered), we require that the m minimization conditions, 


/l = /2 = ■ ■ ■ = /m = 0, 


( 27 ) 


continue to hold under an arbitrary variation in a model parameter p —)• p + (5p, so that the 
variations 6 fi satisfy 

Sfi = Sh = --- = Sfm = 0. (28) 

Each fi is a function of the VEVs Vj and I running parameters evaluated at the scale of 
EWSB, fi = fi{vj, Qk). Thus for each fi we hnd that 


X ^ ^Ji X ^ _ g 

dvj dp dqk dp 


(29) 


The quantities ^ are the elements of the CP-even Higgs squared mass matrix of the 
model before rotating into the mass eigenstate basis. When evaluated for all n model 
parameters, the above system of equations can be concisely expressed as 


Ml 


^771 

dvi ' 



. ^ \ 



. \ 

dpi 

dp„ 


dqi 

dqi 


dpi 

dpn 

dVm 

dVm 


dfm 

dfm 


dqi 

dqi 


(30) 


The quantities forming the hrst matrix on the right-hand side, along with M|, are eas¬ 


ily calculated by differentiating the conditions in Eq. (16) with respect to the VEVs and 


the running parameters. The remaining derivatives dqk/dp must be determined using the 
RGEs. Once these have been obtained, it is straightforward to solve for the dvi/dp. The 
sensitivities Ap are then simply linear combinations of the dvi/dp and dqk/dp. The effects of 
radiative corrections may be easily included by including the Coleman-Weinberg potential 
contributions AV in the EWSB conditions. Here we use the one-loop corrections given in 


Eq. (14) 


Evaluating the derivatives dqk/dp must, in general, be done by numerically integrating 
the two-loop RGEs. This is time consuming and presents an obstacle to doing large scans 
of the parameter space. For studying models dehned at low energies, as we do here, we can 
take advantage of the fact that the running is over much smaller scales than when evolving 
up to the GUT scale. This makes it possible to use approximate analytic solutions to the 
RGEs that exhibit good accuracy over the range of scales considered. Given the two-loop 
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RG equation for a parameter q, 


^ = B =- 

dt ^ 


A''’ + 




t = In 


Q 


(31) 


( 167 ^ 2 ) 2-9 ’ ^ - Mx 

a Taylor series expansion of the solution may be used to obtain the parameter at the scale 


(liQ) = q{Mx) + / Pqit') dt' ^ q{Mx) + 






( 2 ) 


( 1 ) 




+ 


d/?, 

327r2 dt 


+ 0{t^). (32) 


Expanded to this order, we obtain the leading log (LL) and next-to-LL (NLL) contributions 
at two-loop order. The 0{t^) terms not displayed above are formally of three-loop order and 
are neglected. The derivative of the one-loop (3 function is given by 


dA 


( 1 ) 




dt Ibvr^ ^ dqk ’ 

Qk 


where the sum is over all running parameters appearing in /3q^\ The A functions 


(33) 


appearing 


on the right-hand side of Eqs. (32) and (33) are evaluated at the scale Mx, giving a simple 
analytic expression for the parameters at the scale of EWSB in terms of the model parameters 
at Mx- Explicit results for the relevant series expansions in the MSSM and Eg models are 
presented in Appendix 


V. RESULTS 

Using the approach outlined above, we are able to scan the low-energy parameter space 
of the MSSM and EgSSM and calculate the hne-tuning in each. To do so, we implemented 
the above expressions for computing the hne-tuning in a modihed version of the EgSSM 
spectrum generator that was used in Ref. [I6]. This code implemented two-loop RGEs for 
all parameters except the soft scalar masses. In order to properly include the hne-tuning 
impact of the SU{3) gaugino soft mass M 3 , we have extended the original code to make use 
of the two-loop RGEs generated by SARAH |140H143j and FlexibleSUSY jl44j . which also 
makes use of SOFTSUSY |145[ 1146] . The GP-even Higgs masses are calculated including 
the leading one-loop ehective potential contributions given in Ref. [H] and for the light 
Higgs we use the leading two- loopj^ contributions from Ref. [13] which are a generalization of 

^Two-loop corrections calculated for a nonminimal SUSY model may now also be obtained from SARAH 
[14711148j . However, this was not available when the numerical work for this paper was carried out, and 
such corrections go beyond the required precision for studying fine-tuning here. 
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the corrections in the MSSM and NMSSM calculated using effective field theory techniques 
[M USD]. To scan over the MSSM parameter space, the equivalent MSSM fine-tuning 
expressions were implemented into a modified version of SOFTSUSY 3.3.10 [145] . For 
consistency with the results produced in the Eq models, and for computational speed, for 
our main scans only the dominant one- and two-loop corrections to the CP-even Higgs 
masses were included. Finally, in all of the results below the fine-tuning was evaluated at 
the scale Q = Msusy = where mi^ ^ are the running DR stop masses evaluated at 

Q = Msusy- 

As discussed in the Introduction, many recent papers interested in natural supersymmetry 
have focused on light stops, with much theoretical effort to find models where it is easier 
to get a 125 GeV Higgs boson and light stops simultaneously and much experimental effort 
to search for light stops. This is entirely appropriate since there are many good reasons to 
expect the soft masses to be set at high energies. However, that is not the only possibility 
and the fine-tuning problem depends strongly on the RG evolution from the GUT scale, as 
the soft Higgs masses that appear in the EWSB conditions pick up contributions from the 
soft squark masses. 




FIG. 1. Left panel: Scatter plot of fine-tuning in the MSSM as a function of the lightest stop mass, 
for the cutoff scales (from bottom to top) Mx = 20 TeV, Mx = 50 TeV, Mx = 100 TeV 
and Mx = 10^® GeV. Right panel: Scatter plot of hne-tuning in the MSSM as a function of the 
lightest Higgs mass, nih ^, for the cutoff scales (from bottom to top) Mx = 20 TeV, Mx = 50 TeV, 
Mx = 100 TeV and Mx = 10^® GeV. 
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To illustrate this, in the left panel of Fig. we show the variation in fine-tuning for 
Mx = 20 TeV, 50 TeV, 100 TeV and 10^® GeV when we scan over the stop masses and 
mixing, with 500 GeV < < 10 TeV and —3810 GeV < At < —20 GeV. The 

remaining parameters we fix such that at Msusy they have the values /i = —97.5, B = —84.8, 
Ml = 92.1, Ms = 95.9, M 3 = 352, A/, = -117.9, A^ = -7.8, tul, = 400, rUe, = 204, 
mQi 2 = 438, mut ^2 = 436 and = 438 GeV (i = 1,2,3). Here we denote by Mi, Ms 
and M 3 the soft gaugino masses for 1/(1), SU(2) and SU(3), while A), and A^- are soft 
trilinear couplings and the are soft scalar masses for the indicated fields. The soft 
bilinear coupling B is dehned such that at tree level the mass of the GP-odd MSSM Higgs 
boson reads m\ = 2B^/ sin 2 /?. All off-diagonal couplings and scalar masses are set to zero, 
as are the first- and second-generation Yukawa couplings and soft trilinears. Although we 
should stress that making this choice will lead to a spectrum which is in conflict with the 
LHG limits, doing so ensures that fine-tuning due to the other parameters is small, so that 
we avoid washing out the hne-tuning impact of the stops when the tuning is smalQ as can 
be the case when the stop masses are less than 1 TeV. Note that the Higgs mass is also 
allowed to vary in this scan, as shown in the right panel of Fig. [Tj This illustrates the tuning 
problem which people have been worrying about since the discovery of the 125 GeV Higgs 
boson as we see that raising the stop masses is also pushing up the Higgs mass, meaning 
that heavier Higgs masses require more hne-tuning. However, for a low value of the UV 
scale this tuning is not so severe unless the stops are very heavy, and a 125 GeV Higgs can 
be obtained without much tuning in this unrealistic case where we have minimized other 
sources of tuning. On the other hand, the tuning becomes more severe as we increase the 
cutoff such that for Mx = 10^® GeV a lightest stop mass of 1-3 TeV can result in a hne- 
tuning of ~ 100 - 1000 and the minimum tuning we hncQfor a 125 GeV Higgs is ~ 200, as 
shown in Fig. 

Since the stop mass does not have such a large impact on the hne-tuning when the cutoh 
scale is very low we can use this to see more clearly the impact of the Z' mass on hne- 
tuning. To do so we select a hxed low cutoh of Mx = 20 TeV and compare the hne tuning 
between the MSSM and EgSSM for two diherent values of the Z' mass. We choose to look 
at Mz' = 2.5 TeV, which is just above the current limits, and Mz' = 4.5 TeV, which should 

^For models in which the spectrum is heavier, when the stop masses are small the fine tuning reaches a lower 
bound imposed by other heavier parameters. 

^Note that in the calculation of the Higgs mass there is a significant theoretical error, even with leading two- 
loop corrections, which should be considered when thinking about what the results imply for the minimum 
fine tuning in the model consistent with the recent]^scovery of a 125 GeV Higgs. 



be in reach in run II at the LHC |151j and then compare the hne tuning calculated in each 
case to the tuning in the MSSM. For this, we have performed a six-dimensional parameter 
space scan in both the MSSM and EgSSM, varying those parameters most relevant for the 
hne tuning and the Higgs mass. Therefore, the set of parameters which we vary includes /i, 
B and tan/3 for the MSSM, and A, A\ and tan/3, for the EgSSM, which appear at tree level 
in the EWSB conditions of the models. While the RGE contribution from large stop masses 
to the hne tuning is small for such a low cutoh scale, the stop contributions to the ehective 
potential can play a signihcant role in reducing the hne tuning. For this reason it is still 
important to properly treat the tuning associated with stop contributions to the one-loop 
ehective potential, and so we also scan over the soft masses and the stop mixing 

A-t. The relevant parameters and ranges that were scanned over are summarized in Table |TT| 
In addition to this we also repeat each scan for three diherent values of M 2 to allow more 
variation in the chargino masses. 

In this case, we now consider realistic scenarios, where the parameters that are not 
scanned over are set to values which keep the associated states comfortably above their 
experimental limits. So in both the MSSM and EgSSM, all other soft scalar masses are set 
to 5 TeV. We require a valid spectrum with no tachyonic states to exclude points which 
would have an unrealistic minimum, for example due to the appearance of charge or color 
breaking (CCB) minima. We work in the third family approximation, taking the hrst- and 
second-generation Yukawa couplings to be zero, and we also assume that their associated 
soft trilinears vanish. Similarly, we take Ai, = At- = 0 GeV. The U{1) gaugino soft mass Mi 
was hxed to Mi = 300 GeV, and we hx M 3 = 2000 GeV. Additionally, in the EgSSM the 
U{1)n gaugino soft mass M[ is held hxed with M( = Mi = 300 GeV, and /i' = 5 TeV. 

In Fig. results from the scan are plotted showing the tuning for each case against the 
lightest Higgs mass. As expected, the dependence on the Higgs mass is now quite weak, 
while the minimum tuning in the model for the EgSSM is increased by the mass of the Z' 
boson. So in the case of a very low cutoff the tuning required to get a 125 GeV Higgs is 
not so large. However, the tuning from the Z' mass appears already at tree level and is, 
therefore, not suppressed when the cutoff scale is low. In our scan we hnd that, for the 
points satisfying the current limit on the mass of the Z' boson and having an approximately 
125 GeV Higgs, the minimum hne-tuning that can be achieved is Amin ~ 121. If run H of 
the LHG further pushes up the limit on the Z' mass to be above 4.5 TeV then the hne-tuning 
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MSSM 

EeSSM 

2 < tan d < 50 

2 < tan /? < 50 

— 1 TeV < /X < 1 TeV 

-3 < A < 3 

-1 TeV < B < 1 TeV 

-10 TeV < Aa < 10 TeV 

200 GeV < mQ 3 < 2000 GeV 

200 GeV < mQ 3 < 2000 GeV 

200 GeV < m „3 < 2000 GeV 

200 GeV < m „3 < 2000 GeV 

-10 TeV < At <10 TeV 

-10 TeV < At <10 TeV 

Ma = 100, 1050, 2000 GeV 

Ma = 100, 1050, 2000 GeV 


TABLE II. The parameters scanned over and the ranges of values used in the MSSM and the 
EeSSM models. 

in the model will be greater than at least Amin ~ 394 for a Higgs mass between 124.5 and 
125.5 GeV. 

This demonstrates two important points about these U{1) extensions-hrst, that limits 
on the Z' mass play an incredibly important role in constraining natural scenarios in such 
models and, second, that the tuning from the Z' limits in these models depends less on 
assumptions about SUSY breaking than the tuning required by the 125 GeV Higgs mea¬ 
surement which concerns people in the MSSM. 

There is another limit which plays a similar role. Ghargino limits directly constrain the 
/X parameter (or effective /i parameter in these U{1) extensions). The LEP bound |152] on 
chargino masses, excluding m-± < 104 GeV, implies that |/x| should only be greater than 
~ 100 GeV, which is not substantially larger than Mz- Gonsequently the bound from LEP is 
not high enough to have an impact on the hne-tuning obtained in the models and parameter 
space regions that we have studied, as we have checked explicitly. Signihcantly larger lower 
bounds on the /x parameter, and therefore on the hne-tuning, may arise from chargino limits 
coming from LHG searches. However, the chargino limits from the LHG depend on whether 
there are light sleptons or sneutrinos and the mass difference between the lightest chargino 
and lightest neutralino. Gurrent limits placed by GMS and ATLAS extend up to m^± 
700-740 GeV if there are light sleptons |153l I154j with much weaker bounds if there are no 
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FIG. 2. Scatter plot of fine-tuning vs lightest Higgs mass for the MSSM (light blue, bottom band), 
EgSSM with Mz' = 2.5 TeV (dark blue, middle band) and EgSSM with Mz' = 4.5 TeV (dark 
yellow, top band). Note that there are points for which the hne-tuning in the MSSM and EgSSM 
with Mz' = 2.5 TeV is larger than is visible on this plot and those below; however, these points 
are obscured by the overlaid data for the EgSSM with Mz' = 4.5 TeV, and it is the lower bound 
on the achievable tuning that is of interest here. 

light sleptons or sneutrinosj^ 

Nonetheless, for the MSSM the impact of potential chargino mass limits is shown in 
Fig. 1^ There we see that if the full parameter space with m-± < 700 GeV was excluded, the 
impact would be to make the tuning in the MSSM with a 20 TeV cutoff similar to that of the 
EgSSM with the same cutoff and a Z' mass just larger than current limits. In the EgSSM, 
while raising the chargino limit can have the same impact in principle, due to current limits 
on the Z' mass already imposing a signihcant degree of tuning, chargino masses do not make 
much of a noticeable change. 

^Useful summary plots of these limit may be found on the public pages of ATLAS, https: 
//atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/CombinedSuimnaryPlots/SUSY/ATLAS_SUSY_EWSummary/ 
ATLAS_SUSY_EWSuinmary.png and CMS http://cms.web.cern.ch/sites/cms.web.cern.ch/files/ 
styles/large/public/field/image/lmage_03_exclusion_Combined.png?itok=8FMBpu_l, 
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FIG. 3. Scatter plot of fine-tuning vs lightest Higgs mass in the MSSM with 200 GeV < m-± < 400 
GeV shown in light blue (bottom band), 500 GeV < m-± < 600 GeV in dark blue (middle band), 
and 700 GeV < m-± < 800 GeV in dark yellow (top band). 


The exact level of tuning from the Z' depends on the charges of the extra U{1) gauge 
symmetry it is associated with. In Fig. |^we look at the fine-tuning for other U{1) extensions 
for the same Z' masses as we did for the EgSSM. To simplify the analysis we fix tan/? = 10, 
but scan over the remaining parameters as in Table and fix the rest to the same values 
we did in the scan carried out for Fig. In order to more clearly identify the lower bound 
on the obtainable tuning in each model, the parameter values for points in these main grid 
scans with a low fine-tuning were then used as the starting points for smaller scans about 
those values. In these smaller scans the parameters were more hnely varied to populate the 
low fine-tuning regions. 

As can be seen in Fig. the severity of the tunings varies quite a bit. This is because 
the charges appear as coefficients in front of the Z' mass in the EWSB condition. These 
charges change the value of the coefficient d in Eq. ( [2T| ). The values of the coefficient d in 
each model, for tan/? = 10, is {—0.01, 0.40, 0.50, 0.81} for {f/(l)/, U{1)n, U{l)n} and 
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FIG. 4. Top left panel: Scatter plot of the fine-tuning vs lightest Higgs mass in the 17(1)/ model. 
Top right panel: Scatter plot of the fine-tuning vs lightest Higgs mass in the U{1)^ model. Bottom 
panel: Scatter plot of the fine-tuning vs lightest Higgs mass in the U (1)^ model. In each plot points 
with Mz' = 2.5 TeV are shown in dark blue (bottom band), and points with Mz' = 4.5 TeV are 
shown in dark yellow (top band). 


this determines which of the models are most tuned. 

Interestingly, the coefficient d is very small (and negative) in the case of the 17(1)/. This 
allows a dramatic reduction in the hne-tuning from the 17(1)/ symmetry. This is a result 
of the Hu charge associated with 17(1)/ vanishing, which means that the D terms to the 
lightest Higgs which is predominantly Hu at large tan/? are suppressed, making it difficult 
to raise the Higgs mass in the same way as happens in the other models and explaining why 
heavier Higgs values in this model can’t be obtained. Therefore, the hne-tuning behaviour 
in this model is closer to that of the MSSM, and in this case raising the Z' mass limit to 4.5 
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TeV will have little impact on naturalness. From naively estimating the tuning, using the d 
coefficient one can estimate that Z' limits need to be around 15 TeV before they will raise 
the tuning in this model. 

Finally we want to emphasize that while in Fig. |^the EgSSM looks more hne-tuned than 
the MSSM this depends on the high scale boundary, Mx, where the parameters are assumed 
to be set by some SUSY breaking mechanism. Indeed in Ref. [16] a constrained version of 
the EgSSM, with the high scale boundary at the GUT scale, is considered and there the 
cEgSSM was found to be less tuned than the cMSSM. Since a 125 GeV Higgs can be achieved 
in the EgSSM with lighter stops, then if the cutoff is large, the larger stop masses of the 
MSSM can make that model more hne-tuned due to large RGE effects. 

To further illustrate this point, we looked at how the tuning varies with Mx for low 


tuning benchmarks in the MSSM and EgSSM. These benchmarks are dehned in Table III 


and the results are shown in Fig. Since the behavior is quite complicated we now discuss 
these in detail as it provides some insight into the many differences in the tuning between 
the two models. 

In the top panel one can see that the MSSM BMl tuning (dotted curve) steadily climbs 
as the cutoff scale is increased, as one would expect when the tuning originates from large 
soft masses entering from the RGEs. The panel on the middle left conhrms this, showing 
that the largest tuning contributions come from and A ^2 with the former being the 
larger sensitivity until Mx ~ 10® GeV at which point /S .^2 takes over, leading to the small 

Hu 

kink in overall tuning which can be seen in the dotted curve in the top panel. In this case 
we have chosen a point with large mixing, which is known to reduce the MSSM tuning. We 
found this does not eliminate the tuning as there is still a strong sensitivity to At, but we 
did hnd that large mixing lead to less hne-tuning overall for the points we examined. 

Gomparing the MSSM tunings to the EgSSM tunings one can see that which point is 
more hne-tuned depends on the scale at which the parameters are dehned. This illustrates 
that any statement about which model is more tuned depends on the high scale boundary, 
Mx. 


For EgSSM BMl the hne-tuning is shown by the solid curve in the top panel of Fig. 
and the individual sensitivities are given in the middle right panel. The tuning actually 
reduces initially as the cutoh is is increased from 20 TeV. This occurs because the largest 
sensitivity is initially A^ (shown in solid light blue in the middle right panel). This contains 
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FIG. 5. Top panel: Scatter plot of the fine-tuning as a function of the cutoff scale Mx for the four 


benchmark points given in Table III, Middle left panel: Individual sensitivities for MSSM BMl 
plotted against the high scale Mx which give the overall tuning shown by the dotted line in the 
top panel. Middle right panel: Individual sensitivities for EgSSM BMl plotted against the high 
scale Mx which give the overall tuning shown by the solid line in the top panel. Bottom left panel: 
Individual sensitivities for MSSM BM2 plotted against the high scale Mx which give the overall 
tuning shown by the dash-dotted line in the top panel. Bottom right panel: Individual sensitivities 
for EgSSM BM2 plotted against the high scale Mx which give the overall tuning shown by the 
dashed line in the top panel. 
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MSSM BMl 

MSSM BM2 

EeSSM BMl 

EeSSM BM2 

tan l3{Mz) 

10 

10 

10 

10 

s(Msvsy) [GeV] 



6700 

6700 

«1,2,3(.^SUSy) 



0.6 

0.52 

'^1.2(AfsuSY) 



0.2 

0.13 

MefflA^SUSY) [GeV] 

689.7 

1013.5 

1093.3 

1313.0 

Beff(MsusY) [GeV] 

345.7 

1032.5 

3792.7 

817.8 

AriMsvSY) [GeV] 

0 

-5057.9 

0 

-88.5 

Ab{MsvSY) [GeV] 

0 

-5707.2 

0 

-1720.7 

At(MsvSY) [GeV] 

-3335.7 

-2734.8 

-1100 

-1103.2 

mii jlA^suSY) [GeV^] 

2.5 X lO’’ 

6.35 X 10® 

2.5 X 10'^ 

4.94 X 10® 

^I^I^SVSy) [GeV^] 

2.5 X lO’’ 

6.22 X 10® 

2.5 X 10'^ 

4.90 X 10® 

^ei^al-^SUSY) [GeV^] 

2.5 X lO’’ 

6.27 X 10® 

2.5 X 10'^ 

5.21 X 10® 

"iesl-^SUSY) [GeV^] 

2.5 X lO’’ 

6.03 X 10® 

2.5 X 10'^ 

5.11 X 10® 

”^Qi, 2 (^SUSy) [GeV^] 

2.5 X 10^ 

7.37 X 10® 

2.5 X 10'^ 

5.76 X 10® 

mQjlA^suSY) [GeV^] 

4.45 X 10® 

3.97 X 10® 

4.50 X 10® 

3.61 X 10® 

™Si, 2 (^SUSy) [GeV^] 

2.5 X lO’’ 

7.30 X 10® 

2.5 X 10'^ 

5.54 X 10® 

m-i^iMsusY) [GeV^] 

4.0 X 10® 

6.60 X 10® 

5.86 X 10® 

2.04 X 10® 

m|i_ 2 (AfsuSY) [GeV^] 

2.5 X lO’’ 

7.30 X 10® 

2.5 X 10'^ 

5.88 X 10® 

■ml^{MsvSY) [GeV^] 

2.5 X lO’’ 

7.03 X 10® 

2.5 X 10'^ 

5.78 X 10® 

^Hji^SVSY) [GeV^] 

1.82 X 10® 

8.96 X 10® 

4.06 X 10'^ 

1.04 X 10^ 

™lf„('^SUSY) [GeV^] 

-3.60 X 10® 

-9.35 X 10® 

5.0 X 10® 

-2.66 X 10® 

”i|(A^susy) [GeV^] 



-3.10 X 10® 

-3.17 X 10® 

Mi{Msusy) [GeV] 

300 

260.8 

300 

173.4 

M 2 (Msusy) [GeV] 

2000 

479.2 

1050 

281.4 

M3(Msvsy) [GeV] 

2000 

1312.3 

2000 

1200 

M[{Msvsy) [GeV] 



300 

175.2 

Mz< [GeV] 



2473.2 

2512.7 

m^i [GeV] 

124.3 

124.4 

125.0 

126.2 

[GeV] 

1942.1 

861.6 

993.8 

1665.0 

[GeV] 

2220.1 

2023.9 

1174.8 

2094.4 

rrig [GeV] 

2259.8 

1472.9 

2290.0 

1407.4 

A{Mx = 20 TeV) 

157.3 

242.8 

165.3 

402.1 

A(Mx = IQi® GeV) 

1089.0 

949.0 

1722.3 

546.7 


TABLE III. Parameters for the MSSM and EgSSM benchmark points. In the EgSSM, we define 
/Teff = As/\/2 and Beg = Ax- The soft masses and are those that satisfy the 

EWSB conditions including one-loop Coleman-Weinberg corrections involving the top and stops. 
For EgSSM BMl (BM2) we also set n' = 5000.0 (897.9) GeV, = 5000.0 (—4.21 x 10^) GeV^, 
^«i. 2.3 = 0 (-1389.2) GeV, = 0 (-52.9) GeV, = 2.5 x 10^ (4.81 x 10®) GeV^, 

= 2.5 X 10^ (4.90 x 10®) GeV^, = 2.5 x lO’^ (4.46 x 10®) GeV^, = 2.5 x 10^ 

(4.81 X 10®) GeV^, = 2.5 x 10^ (5.28 x 10®) GeV^, mjj, = 2.5 x lO'^ (4.94 x 10®) GeV^ and 
= 2.5 X 10^ (4.87 x 10®) GeV^. 
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some terms proportional to M'^, which provide the dominant contribution to this sensitivity 
at very low Mx- However, as Mx is increased contributions from the soft masses become 
more important and these actually start to cancel the large contribution to coming from 
Mz' until Aa passes through zero. At the same time though these large soft masses also 
cause other sensitivities to grow, in particular A^g- The hne-tuning rises with Mx once 
Mx ^ 10^ — 10® GeV, but remains lower than that of the other points, until Mx ~ 10® 
GeV. Eventually the Am 3 sensitivity leads to this point being the most hne-tuned of the 
four shown in Fig. 

Although the gluino mass and M^{Msusy) have similar values to those in the MSSM 
BMl point, in the EgSSM M^{Mx) is larger due to the altered RGE running from exotic 
mattei]^ This is why this EgSSM BMl has a larger tuning at larger values of Mx, coming 
from Am 3 - 

Interestingly other sensitivities are suppressed by this effect since at the same time larger 
M 3 at higher scales reduces the soft squark masses at Mx- Therefore, the stop mass contri¬ 
butions are ameliorated, compared to the MSSM, both by allowing lighter stops at Msusy 
and by the modified RGE running. Nonetheless the stops still do lead to A ^2 increasing 

Hu 

with the cutoff through the usual mechanism^ 

By contrast the tuning for EgSSM BM2 is very different, as is shown by the dashed line in 
the top panel, with the individual sensitivities given in the bottom right panel. This point 
was chosen as it had a much lighter gluino mass that is just above the experimental limit of 
1.4 TeV m- At 20 TeV this benchmark is not amongst the lowest tuned points, since at 
that scale the tree-level tuning from Mz' dominates. However, the reduction in M 3 means 
that Ams is substantially lower and only becomes the dominant tuning at a much larger 
scale of Mx ^ 10^^ — 10^® GeV, giving a tuning at 10^® GeV of ~ 546, which is far below 
that of the other three benchmark points. 

In addition to this, the soft parameters in EgSSM BM2 follow a pattern similar to that 
found in the constrained model. With the exception of the parameters rnq^, 


and M3, the values of which are given in Table HI, the soft masses at the SUSY scale 


correspond to the values that result in the cEgSSM with mo = 2.2 TeV, M 1/2 = 1003 GeV, 

®This altered RGE running is a result of the exotic matter introduced to keep the extra U{1) anomaly free. 

’^Wherein m\^{MsuSY) receives a positive contribution to it’s mass from rrijj^{Mx) and a negative con¬ 
tribution from mQ^{Mx) and mf^^{Mx), allowing heavy stop masses to cause fine-tuning. In this case 
m\^{MsusY) is held hxed so as the cutoff increases the values of these soft masses at Mx will be larger 
and there will be a bigger cancellation between them, increasing the sensitivity of Mz to both and the 
soft scalar masses for the stops. 
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Aq = 500 GeV, = 0.1923, X{Mx) = 0.2646 and \i^ 2 {Mx) = 0.1. This leads to a 

significant reduction in the contributions to the RG running of and ttiq^ coming from 
terms of the form and, to a lesser extent, Here we dehne for the EgSSM (see 


also Eqs. (B8 B9) for general U{1) inspired models) 


d 

f 2 o 2 | 2|2 2,2 2,2 25 2,2 

El = 2^ i^mQ. - 2m„^ + m^. + m^. - + m^, - + m^, 

i=l 

3 

E) = ^ + 3m^, + 6m^. + + 4m|, — 4m^„ — 6m^d + 5m|, — 9m^ — 6m] 


+ 4m^, - 4m^. 


In the unconstrained case, this contribution acts to drive up the values of rug^ and m'jj^, 
and thus the associated tuning sensitivities, at the cutoff scale Mx- In the case of EgSSM 
BM2, on the other hand, the reduced splitting between the soft masses leads to a much 
smaller contribution from these terms. Together with the reduction in M3 described above, 
this allows to maintain the observed low hne-tuning at very large values of Mx- 

MSSM benchmark BM2 (dash-dotted in top panel, individual sensitivities in bottom left 
panel) is designed to be similar to EgSSM BM2, for a reasonable comparison. However, from 
the individual sensitivities one can see that the behavior is quite similar to MSSM BMl, 
though in this case A^2 becomes the largest tuning at a higher Mx and does not reach 

Hu 

such large values, since more of the tuning is from the mixing in this case. 


VI. CONCLUSIONS 

Prior to stringent experimental constraints on the mass of the lightest Higgs boson and 
squarks in supersymmetric models, a simple picture of a natural SUSY model emerged 
from theoretical reasoning, with soft masses set to similar values at the GUT scale through 
local gravity interactions with the hidden sector. Through the use of renormalization group 
running, one can then see that at the EW scale the stops enter the EWSB condition for 
Mz] therefore, it was expected that these masses should not be much bigger than 100 GeV. 
However, to disturb this elegant picture, hrst LEP placed constraints on the Higgs mass, 
requiring it to be above 114.4 GeV |156[ I157j . which already introduced signihcant tuning 
for constrained models since heavy stops are required to raise the lightest Higgs mass above 
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its tree-level bound. Then, recently, this problem got much worse since the LHC measured 
the Higgs mass to be around 125 GeV. 

17(1) extensions motivated by the p problem, Eq GUT theories and the connection to 
string theory contain both F- and D- term contributions to the light Higgs mass which 
can raise the tree-level mass, evading the need for large radiative corrections to increase it. 
However, such models come with their own fine-tuning problem, where the Z' mass appears 
in the EWSB condition for Mz at tree level. While in a previous study of the constrained 
EgSSM it was found that the tuning is less severe than the MSSM, it was still significant. 

In light of such difficulties it is worth considering whether the simple picture which 
emerged is wrong in some way and if there are other possibilities that allow naturalness. Or 
to phrase this in a more challenging manner are there ways to constrain the naturalness of 
these models that do not rely upon assumptions about how SUSY is broken? 

We have investigated this question here in the context of the MSSM and 17(1) extensions. 
Since the RG evolution links the soft masses together and causes these problems from stop 
and gluino masses the most conservative approach to placing naturalness limits is to choose 
a low cutoff. We find that in the MSSM the most direct way to constrain naturalness in the 
model without making assumptions about the SUSY breaking scale is through limits on the 
chargino masses. Gurrent LHG limits on charginos are not model independent and thereby 
leave many gaps where one can have light charginos. 

In contrast we find that in 17(1) extensions of the MSSM there is an additional way to 
constrain the naturalness of the models, which is through the Z' mass limit. We hnd when 
we impose a low cutoff of 20 TeV for setting the soft masses, the lowest tuning in the EgSSM 
compatible with a Z' mass of 2.5 TeV was A 121, while if the LHG run H can place a limit 
of 4.5 TeV on then the tuning would be approximately 394. By comparison the current 
situation in the MSSM only requires a tuning of around 38. This should be interpreted as 
saying that in the most conservative limits one can place on naturalness in these models, 
the tuning in the EgSSM is worse. However, if there are no charginos below 700 GeV then 
the situation in the two models would be the same. 

This should also be contrasted with what happens as we raise the high scale boundary, 
Mx- We showed that for our benchmark points, which one is more tuned depends very 
strongly on Mx- The EgSSM tuning is sufficiently complicated by the interplay of these 
different sources of tension in the EWSB conditions that a small reduction in hne-tuning 
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can even occnr for a moderate increase in Mx- However, as Mx increases towards the scale 
where the gange conplings nnify, the familiar tnnings do dominate, thongh with tnnings from 
the glnino mass appearing to be more signihcant relative to those from soft scalar masses. 

We also looked at the tnning in different U{1) extensions for hxed tan/3 = 10. We fonnd 
that in every case except for the U{l)i the hne-tnning was mnch worse for the larger Z' mass, 
farther emphasizing the importance of this in U{1) extensions. The U{l)i model showed 
the least tnning dne to the vanishing charge of the state. This model is qnite interesting 
in the sense that it provides a solntion to the /i problem while avoiding the large tnning 
(with cnrrent limits) from the Z' mass. However, one shonld remember we are looking at 
conservative limits on natnralness here and there is no solntion to the nsnal tnning coming 
from the large stops needed to get a 125 GeV Higgs in this model, which will be a problem 
as the UV scale is raised. 
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Appendix A: FINE-TUNING MASTER FORMULA 


To write down the tree-level master formula, it is convenient to dehne the quantities 

dfj dfk 


— ^ijk 


ds dtanfS 


with /i,/ 2,/3 as given in Eq. (16). The relevant partial derivatives are 


dfi 


Stan f3 


‘^Mz 2 a f a i ■ a 
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dfi _ 2Mz 
ds g 

dh 2 Mz 
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df2 _ 2Mz 
ds g 

dh 2M| 


\-%(QiQ2-Ql) 


Ayl, 


d tan [3 g"^ 


s (A^ + giQ 2 Qs) sin/3 - cos/3 
cos^ (3 g'lQss {Q 2 — Qi) sin 2/3 — 72 AAa cos 2/3 


dh 2 2A2M| 


4M| 


(Qi cos^ /3 + Q 2 si7 /3) + 


For a running parameter g appearing in the tree-level EWSB conditions, the corresponding 
contribution to the individual sensitivity can then be written 


X _ dfi , dh dh 

“Z-1“ ^2 ZX-1“ ^3 “Z” • 


(A2) 


dq dq dq 

It is straightforward to compute the appropriate derivatives directly from the EWSB con¬ 


ditions, Eq. (16). Similarly, the quantity C appearing in Eq. (26) is given by 

dh 


C* = z Ui 
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Appendix B: RGE CONTRIBUTIONS 


Provided that one does not run over too large a range of scales, the solutions to the RG 


equations for a model can be reasonably well approximated by a Taylor series, Eq. (32). For 
a parameter p, this reads 


q{Q) = q{Mx) 


t 


Ibvr^ 


4" + 


( 2 ) 


Ibvr^ 


+ 




(Ibvr' 


2\2 Q 


bf\M_ 


X} 


where we have for convenience dehned 

Mx 

We have constructed the necessary series solutions in both the MSSM and the U (1) extended 
models. Due to the smallness of the Erst- and second-generation Yukawa couplings, we 
neglect them in our calculations. The corresponding soft SUSY breaking trilinears are 
likewise omitted. Additionally, all soft mass matrices are assumed diagonal, and the gaugino 
masses are taken to be real. 

In the MSSM, the relevant parameters for the hue-tuning calculation are /x, B, 
at tree level. For the renormalization group running of the relevant parameters SOFTSUSY 
uses the one- and two-loop RGEs from |15811159j . The corresponding 0{t^) contributions 
are 


ft(2) = ^ 
^ 2 
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In these expressions the qnantity S is dehned by 

3 

= ^Hu -^Hd + Y^ +^l)- (B2) 

i=l 

If, in addition, the one-loop contributions to the effective potential from top and stop loops 
are included, it is also necessary to construct the expansions for rriq^, m\^ and At. The 
coefficients read 


= 2 %f + ^la + 2 A^) + 247/^^ {mjj^ + m%^ + 771^3 + 2Al) 
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+ 96^3^M| - —<7^ (5 + AMl) , (B3b) 

~ + 24?/^y4b + lAylyl{At + A;,) + 2ylyl{Ab + A,-) 
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- 64 ^^|/t 2 (Ai - M3) - 3Qgly^{At - M2) - —gly^{At - Mi) 

32 14 

- y^ 73 \^(A - M3) - 67722/2(A - M2) - -gjyliAb - M,) 

572 

- 64gtM, + 12gtM2 + ^g^M^ (B3c) 


(2) 

We can similarly obtain the two-loop f3 functions and coefficients bp ^ for a general set of 
U (1)' charges. Two-loop RGBs for the gauge and Yukawa couplings, gaugino masses and soft 
trilinears, along with the one-loop RGBs for the soft scalar masses, were originally obtained 
in Ref. [IH] for the particular case of the BgSSM. BlexibleSUSY uses full one- and two-loop 
RGBs from SARAH, which for the models considered here are based on Ref. usa and the 
recent extensiorj^ in Ref. |16nj to include models with multiple G(l) gauge groups, in the 
most general case where the trace of the matrix formed from the charges Qj of the U{1)y 
gauge symmetry and Qi of the extra U{1)' symmetry does not vanish, i.e. Yhi QjQi 7^ 0- 
When this trace is nonzero it will also induce gauge kinetic mixing to be generated during 
RGB evolution and this is the case in the models we consider herej^ However, when these 
models are evolved down from the GUT scale, the radiatively generated gauge kinetic mixing 


®In the version of SARAH which we used the extra terms from this extension were included for all terms 
except the trilinear and bilinear soft masses. We have been in contact with the SARAH author about this 
and understand they will be included in future versions. 

®In Eq inspired models with only complete 27-plet matter multiplets this trace would vanish. However since 
we assume some incomplete multiplets so that our models are consistent with gauge coupling unification 
this trace doesn’t vanish. 
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gives an off-diagonal gauge coupling, gn, of just ~ 0.02 at the electroweak scale [13] and so 
it does not play a large role. At the same time if gauge kinetic mixing is included the RGE 
expressions become very large and unmanageable, so we neglected the gauge kinetic mixing 
by setting the SARAH flag NoUlMixing to true. 

At tree level in the EWSB conditions the parameters that must be considered are A, A\, 
'^‘s S'!) 5'2 g[. Neglecting kinetic mixing, the two-loop (3 functions for 

the relevant gauge couplings read 

= f st (B4a) 

= ^9? (SOsfnJ + llTgJ + 135g= + 300g= - lOE, - 15Ej - 65g? - 35g? - 45g3 , 

(B4b) 
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- IOEa (Ql + Q? + Ql) - 30yl (Ql + Ql + Ql)). (B4f) 


In order to keep these expressions compact, we have used the notation 
Sq = Q? = ^ cos^ d+‘^ sin^ 9 + sin 29 
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to denote the trace over the U{iy charges, along with^^ 
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ye 


^q = Y1 = -I^cos0 - -^sm9, 
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nQ = 5^Q 
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2049 


in 483 3^ . ^ 681 2 n ■ 2 n 9 ^ . 3^ 1297 

-cos 9 -I- 1 = cos 0 sm 0 -I-cos 9 sm 9 H- = cos 9 sm 9 -\ -— 

1600 80^15 160 leyiS 576 


siy 9, 


^°The first of these is the trace which is assumed to vanish in Ref. |159j . Although we use the NoUlMixing 
flag to neglect gauge kinetic mixing, SARAH does this by removing the RGE for the off diagonal gauge 
couplings and effectively setting them to zero at all scales by removing all terms involving them from the 
RGEs. Therefore, some terms with this trace remain and the RGEs shown here do not reduce to those 
which one would obtain from Ref. |159j or Ref. [TH| unless Eg = 0. Note, however, these contributions do 
not appear in the corresponding trilinears due to the version of SARAH used. 
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Note that in these expressions the U{1)y and charges are assumed to be GUT- 

normalized. The expressions in terms of the Eq mixing angle 6 follow from the charge 
assignments given in Table and hold provided that G(l) mixing is neglected. Similarly, we 
write 
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It is sufficient for our purposes to consider the trilinear coupling ax = AAa, rather than Ax, 
for which the relevant expressions read 
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- ^g'l {Ql + Ql + Ql) (Sa, - m(Sa) - i 2 ( 7 f (g| + + gy - M{sg 

576 

+ 485 ^ 2^72 +Ml + (g^ + g 2 + g|) , (BTc) 

where 


Sa;^ — AiOAi + \ 2 a \2 + AsOAg, — Kia^i + K2aK2 + ^30k3 5 

Ila^ = AiOAi + X^ax^ + A^OAg, + ^3®K3) 

Note that oa^ = A^Aa-, a* = |/tA, Ob = l/6-4f, and Dehning 

3 

f 2 o2|2|2 2,2 2,2 2A 2,2 

El = 2^ - 2m„, + m^. + + m^u - j - 771 ^, + mj^r, 


S'l — (QQq'm-Qi + + ‘iQdXn'^di + Qe'oil. + 2QLm\_ + 2Q2m‘^u + 2Qim^^d 


+Qsms^ + 3g^^^ + SQnmjj^j + 2QH'mjj, + 2Qjpm^, 
the one- and two-loop (3 functions and the O(t^) coefficients for are 
= 2A2 + m]j^ + m|) + 2a^ + Qyl {mjj^ + m|3 + ml^) + 6al 

Hd 

+ ‘^Vr i^Hd + + 2ar - S^sMl - ^^iMi^ - sgi^f Mf 


-gli:^ + 2Q,g[^j:[, 


(BlOa) 


PU = -S6yt {m]j^ + 772^3 + mX) - 12y* + m^) - 72ylal - 2Aylal 

- ^vWb {rnXd + ^Hu + 2"iQ3 + mX + mX) - 6 {ytat + yhatf 

- 4A^ {mX^ + mX^ + m|) - SA^a^ - [4A- (mX^ + mX^ + 2m| + mXf 

i=l 

3 

+mXu^ + Qk] (mX^ + + 2m| + [4 (A^oa + AoaJ^ 
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+6 {Kittx + {m]j^ + + rrig^ + mlj - 6 {Xat + yta\f 

+ ^‘^alvl + ^Qs + + 2^3) + 32^^ {al - 2yi,abM^) 

+ ^aWt ( 3 ^i„ + ^q3 - 4771^3) - ‘^glyl {llm]j^ - ml^ - Am^ + AM^) 

- tdi («fe - ‘^ybabMi) + \glyl + ml^ + Aml^ + AM^) + X^gl {al - 2y^arMi 


6 


Y1 + 12c/f l/b (Qq + Qd - Q?) 

^=1 


^ + ^Qs + ^X + 2Mf) + 125 (f (Qq + QI — QI) {o-b ~ ‘^yba-bM[^ 

— 2AQig'^yl + Qq^Hq^ + QdXnlX) “*“ ^9iyl {Qt A- Ql — Qi) 

^ i'^Hd A~ ^La + ^63 + 2M(^) + Ag'^ {QI + QI — Ql) {al — 2y^arM[^ 

- SQig^yl {Qimlj^ + QLml^ + Qeml^) - 2AQig'^yl + Quml^) 

+ 4(7f {Ql + Ql — Ql) + ^1 + 2M(^) + 45ff (Q 2 + Q| — Qi) 

3 

X {al - 2XaxM[) - AQig'^ ^ |^2A- (^Qim^d + Q 2 mlju + Qsmf'^ 


i=l 


+3K' (m| + m^^ +m^ 


16 


2 2 (2 O 2 I 2 I 2 

%9i 2^ - 2m„^ + 

2=1 

+ 32Qi5f35ff (2QQmQ^ + Quml. + Qd^nl. + + 

2=1 

3 

29M| + ^ (sm^, + m|, + m^d + m' 

2=1 

3 

2 {Ml + Ml M 2 + M 2 ) + mlj, - ^ + m]ju - ml^d 


+ 3(72 


,2 


' 2 2 
:^2^l 


+ ^2,Qiglgf 
3 


2Qi {M^ + M 1 M 2 + M 2 ) + 4“ 

1 


2=1 
4 2/2 

- -^9,9, 


25^ 


+ 7^9t 891M| + 18m^, 


(sQquiI^ + QimX + Qim\d + Q2r«H^^ 

2=1 
3 

+ 2 ^ (m|^ + - 9m| + 9m^d + 9m|^ + + 28m|^) 

6 Qi( 3 Qci + 3 Q^ — “AQd + 3 Qe — 4 Qi + 3Q2 + Q/77 — Qh' — “AQl 

+ 3Qq — QQu) {Mi + MiMi + Mi^) + 3Q|^?7 i|^ — 3Q|-/?ti^/ 
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4 


+ 3 ^ {Qlm\ + - Qlml. + Q^ml - QWiif + QWi,, - Qlml. 

i=l 

+ Qg^Qi - + ^QiQw^w - ^QiQn’mjj, + Qi ^ Uq^ml. + AQ^m- 

i=l ^ 

- SQDml). - 9Qim]jd + 3Q2mjjn - 9QLm\^ + ^Qqttiq. - 20(5«m^J 

+ SQiQi 3QiMi^ (9Qd 4 4 9 Qd 4 3Q^ + 8Qi + QQ 2 4 4 2Q^/ 4 QQ 

3 

4 ISQq 4 3(5| 4 9(5^) 4 2Q^m^ + 2Q^^,rn\j, + {^Qd'^% 4 3Q^m^_ + 3Q^j^m 

i=l 

4 4 2Q\m\d 4 2Qlm]ju + 2Qim^^ + 6g|m^, + g|m|^ + 

3 

+ 2QiQ‘^m^ + 2QiQ‘]j,m\j, + Qi (3Q^m^^ + 3Q^m^^ + 3g|,m|), + Qlml. 


}l 


i=l 


+ 2Qlm\d 4 2Qlm]ju + 2Qlml^ + 6gJmJ, + Qlm% + 3Qlml^ 


(BlOb) 


b 


( 2 ) 

m2 


= 72|/b^ + 771^3 + 77723 ) + 1447/2a2 + 16//^ ( 777 ^^ + ml^ + 777 ^ + 32ylal 

4 8A^ ( 777 ^^ + 777 ^^ + 777|) + 16A2a^ + Qylyl {m]j^ + 777 ^^ + 2777 J 3 + ml^ + 777 ^ 3 ) 

+ 6 (T/tOf, + yi,atf + 127/2y2 (2777^^ + mg^ + 777^3 + 77713 + ^li) + {Vbar 4 7/^0^)^ 
+ 6A27/2 (r77^^ + 2m]j^ + 777I + 777 q3 + 77723) + 6 (Acj + ytaxf 
4 12A27/2 (2777^^ + m]j^ +ml + mg^ + 777^3) + 12 (Aa^ + T/bOA)^ 

+ AX^yl {2m]jd + ^Hu +^l + 4 777^ + 4 (Aa^ + t/^oa)^ 

3 

+ 4 ^ |^A2A2 ^777^d + mlju + 777^^ + 777^^ + 2777| j + (Aoa, 4 Xittxf 

i=l 

3 

+ ® 4 + 7774J + {\a^. + KiOA)^ 

- 32^17/^2 (777^^ + 777 ^3 + 777^3 + 2M32) - 32^2 (^^2 _ 2 y ^, abM ^) 

- ^^glvl 4 777q3 + 777^3 + 2Ml) - I8gl (al - 2ybabM2) 

- Qghl + 7772 3 + 77723 + 2Ml) - Qgl {al - 27/^a^M2) 

- 6fi'2A2 (777^^ + 777^^ + 777| + 2M22) - 6gl {al - 2\axM2) 

- ^glvl + 2^9 - yl/? (flfe - 27/babMi) 

- {^hd 4 "^L3 + + 2^1) - (flr - 2yrarMi) 
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— + m]j^ + m| + 2Ml) - ^gl (a^ - 2\axMi) 

— ^‘^QiUb {Qi + Qq + Q^d) {^'ud + ^Qa + + 2Mf) — 12g'^ (Q? + Q'q + 

X (a^ — 2yi,abM[') + Qg'iyl {Qi + Qq + Qd) + ‘^Qi'^Qs + 

+ 12Qi5'f(Qi + Qq + Qd) — ^diUr {Qi + Ql + Qe) + ^ea + 2M(^) 

~ ^9l {Qi + Ql + Qe) (®r ~ ‘^yrdrM'i') + 2g^y‘^ (Ql + Qi + Qe) 

X (2Qimjj^ + 2Qim\^ + 2Qim‘l^ + + AQig'^a^ (Qi + Qi + Qe) 

+ 12Qi(y'f|/i (Q 2 + Qq + Qu) + ^Qa + ^na) l^Qigffa^ (Q 2 + Qq + Q«) 

— 45(f (Q^ + Q 2 + QI) + ^5 + 2Mf) — Ag'^ (Q^ + Q 2 + Q|) 

X (a^ — 2\a\M'^ + 2g'^X^ (Qi + Q 2 + Qs) + 4Qi5'f (Qi + Q 2 + Qs) 

3 

X I^A^ 4* + ^5 j + ®Ai + SQifi'f (Qs + Qd + Qd) 

3 

X ^ j + Ok, ~ SSQigfggffMg (2 Qq + Qn + Qd + Qd + Qd) 

2=1 

— 72g2Ml — 12Qig2g'i M 2 (9Qq + 3Qi + 3Qi + 3 Q 2 + Q|^ + Qd') 

— '^9i (^1 + 3M^) — ~^9i9i 4QiM^ {2Qd + 2 Q 2 J + 2Qi) + 6Qe + 3Qi + 3 Q 2 
+ Qjp + Qd' + 3Q i + Qq + 8Q^t) — AM'^ (SQ^ + 3Q^ — 3Q^ + 3Qg — 3Q^ + 3 Q 2 
+ Qw- Qw - 3Qi + 3QJ - 6Qy + (s; - 2 Q 1 S 1 ) S^] - 4Qi(7;" [ 2 Mf (9Q3 

+ 9Q^ + 9Q^ + 3Qe + QQi + 6 Q 2 + 2QM -\- 2Q^, + 6 Qi + 18 Qq + 3Q|' + 9Q^) 

+ (6QiMf - S'l) Sq] . (BlOc) 

Similarly, those for read 


Pill = 2X^ {mjj^ + mjj^ + ml) + 2al + 6yl {m]j^ + + ml^) + 

Hu 

- QglMi - ^-glMl - 8Qlg[^M[^ + ^112, + 2Q^g'^i:\, (Blla) 

+ "^Qa + - ^ytvl + ^Hd + 2^Q3 + mis + 

- 72|/t a? - 6 {ytab + ybatf - 4A^ + m|) - 8A^a^ 

3 r 

- A^ ^ AX] (m]j^ + m]j^ + 2m| + m^d + 

2=1 ^ 

1 ^ 

+ 6k- (m]i^ + m]j^ + 2m| + [4 (A^oa + AoaJ^ 
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+ 6 {niax + ] - 6\%l + m]j^ + m| + rri^Q^ + mX) 

- 2X^yl + ml + m\^ + - 6 (Aa^ + Vhaxf - 2 (Aa^ + 

+ ^‘^dlvt {'^‘hu + ^Qa + ^«3 + 2M|) + 325f| (a^ — 2ytatM^ 

+ + m^3 + + SM^) + (^^2 _ 2ytatMi) 

+ (3m^^ - 771^3 - 2 my + ^-glyl {m]j^ + 771^3 - 2777^ 

+ [^Hf - (""a - 

2 = 1 

+ ^‘^dlUt {Qq + Qu~ QI) {'^Hu + ^Qa + ^U3 + 2Mf) 

+ 12^f (Qj + - Ql) {a^ - 2ytatM[) - 24Q2^fl/t (<52 ^h„ + <5Q"iQ3 + Quml^) 

- 2-^Q^ghb + Qq^Q3 + Qdm%) - ^Q2g'hl + <5L"ii3 + Qe^y 

+ 4(7f A^ (Ql — Q 2 + Q|) + ^77„ + + 2My 

+ 45 'f (Ql — Q 2 + Qs) (oa ~ 2Xa\M[^ 


I + Qoml)^ + Qjfnl^ 


6 

- 4Q25'f ^ |^2A- {Qim\d + Q2mjju + Qs^rilj + 3 k: 

2=1 

+ y (73^1 - 2^n, + dnl + 

2=1 

3 

+ 32Q2^3^f ^ (2 Qq^q, + Q«^^, + Qd^rf, + Q^% + Qoml^^ 

2=1 

3 

29M| + 777^, + 777|y + ^ (^3777 q^ + 777|^ + 777^ + 777^^ 

2=1 

3 

2 (Ml + M 1 M 2 + M 2 ) + m^p - rriff, + ^ ^777|^ - 777 ^^ + mjjn - m^d 


+ 3(72 


2 2 
^2^1 


+ 12Q2(72fi'f 

3 

+ (3<5(3"iQ, 

i=l 
3 

+ ^ + 2^^Di + 10777^ + 54777g, + 18777^ + 4777 q, — 8777^ 

i=l 

4 


2Q2 (-Ml + M^M2 + -M2) + Qii'm^ji + Qjpm 

1 


2 

77' 




891Mi' + I 8777 I 7 


5 


2 /2 


6Q2(3Qi7 + Qtj — 3 Q 15 + 3Qe — 3Qi + 4 Q 2 + Qjjt — Qh' ~ ^Ql 


+ 3Qq — 6Qu) (M 3 + MiM( + M(^) + SQ^pm^p — 3Q^/777^/ 
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4 


+ 3 ^ (Qdm% + - Qlml^ + Qlm\ - Q\m]ja + Qlm]ju - Qlml^ 

i=\ 

3 

+QQ^Qi - + ^Q2 Qw^w - ^Q2QH'm]j, + <52 XI 

i=\ 

-AQoml). + 12Qeml. - SQim^d + 9Q2mjjn - 3QLml. + IQqtjiq. - 4:Quml^ 

+ 8 Q 2 gi 3Q2M^ (9Q‘d + + 3Q^ + QQ\ + 8 Q 2 + 2 Q^ + 2 Q^, + GQj 

3 

+ ISQq + 3Q^g + 9Ql) + 2Q^m^ + 2Q^^,rn\ji + X] {^Qd'^% + 

i=l 

+ 8Qlm% + Qlm^ + 2Qlm%a + 2Qlm]ju + + 6QjmJ, + + 3Q^m^J 

3 

+ 2Q2Q^m^ + 2Q2Q\{im‘]j, + Q 2 X] {^Qd'^di + 


?L 


2=1 


+ 2Qlm]j^ + 2Qlmjju + 2Qlni% + + 3g2rriy 


(Bllb) 


+ ^Q3 + "^na) + 144|/2a2 ^ g^4 ^ ^ ^ ^g^2^2 

Hu 

+ 6i/ii/b + 2mQ3 + ml^ + + 6 {ytab + ybatf 

+ 12A2|/i^ {m]j^ + 2m|^ + m| + m ^3 + 771 ^ 3 ) + 12 (Aa* + T/tOA)^ 

+ 6A^7/2 {2m]j^ + + mlj + 6 (Aa?, + ybttxf 

+ 2A^7/^ ( 2 m|^ + + m| + 777^3 + ml^) + 2 (Aa^ + yrttxf 

3 

+ 4 X + 2777 !j + {\aXi + AjCA)^ 

j=i 

3 

+ ® + "^1,) + 

2=1 

- 32c/2|/2 + 2M|) - 32^2 (^^2 _ 2y^atM3) 

- ^^ght + rriQ^ + mis + 2 ^ 2 ) - 18^2 (a? “ ‘^ytatM2) 

- 6^2A^ + TTil + 2Ml) - Qgl {al - 2\axM2) 

- ^glvl {m-H^ + 777I3 + ml^ + 2 M^) - ^g^ (a? - 2 ytatMi) 

- ^ 1 /?A 2 ( 777 ^^ + 777 ^^ + 777| + 2M^) - ^gl {al - 2\axMi) 

~ ^‘^giUt {QI + Qq + QI) (^h„ + ^Qs + + 2Mf) — 125 (f (^2 + Qq + Qu) 

X {al — 2ytatM'^ + {Q 2 + Qq + Qu) (2^2^27^^ + 2Q2mQ^ + 2Q2ml^ + E'J 
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+ 12Q25'fa? (Q 2 + Qq + Qu) + 12Q2fi'fl/fe (Ql + Qq + Qd) + ^Qa + ^da) 

+ 12Q2fi'fflb (Ql + Qq + Qd) + 4Q2fi'fl/T (Ql + Qi + Qe) + ^ia + ^ea) 

+ 4Q2fl'f flr (Ql + Qi + Qe) ~ ^9l (Qi + Q 2 + Qs) + 2M(^) 

~ ^g'l^ (Qi + Q 2 + Qs) (®A ~ 2\a\M[^ + 25 ff (Qi + Q 2 + Qs) 

3 

+ 4 :Q2g'i {Qi + Q 2 + Qs) ^2 + ®Ai 

2 = 1 
3 

+ 6Q2fi'f (Q 5 + Qd + Qd) ^2 + ^d* + ^Si) 

2=1 

— 96Q2glg'i ( 2 Qq + Qu + Qd + Qd + Qu) ~ ^‘^ 92^2 

— 12Q29291 M 2 (9Qq + 3Ql + 3Qi + 3 Q 2 + Qjp + Qh') + ~^9t (^1 ~ 3M^) 

~ ~^9i9i 4Q2M^(2Q(i + 2Q;^ + 2Q£) + 6Qe + 3Qi + 3 Q 2 + Qjjt + Qu + 3Qi 
+ Qq + 8Qtj) + 4iVf]^^(3Q^ + 3Q^ — 3Q^ + 3Qg — 3Q^ + 3 Q 2 + Q^ ~ Qd' 

— 3Qi + 3QJ - 6Q2) - ( 2 Q 2 S 1 + s;) S^] - 4Q2(7;" [2Mf (QQ^ + 9Q% + 9Q1, 

+ 3Qg + QQf + 6 Q 2 + 2 Q^ + 2 Q^/ + 6 Q|^ + 18Qq + 3Q^ + 9Q^) 

+ (6Q2Mf-S;)Sg], (Bile) 

while those for m| are 


+ ^Hf + "^ 5 ) + 4aA, + 


- SQlg'^M[^ + 2Q5(7;"E;, 


(B12a) 


3 r 

i=l f 

- 32 A 2 a^^ - - I2\^yl (2m^^ + + m| + 771^3 + 771^3) 


- 12\^yl {m]j^ + 2m]j^ + ml + + m%) 


- AX^yl [mjj^ + 2m]j^ + ml + m\^ + ml^) - 12 {Xat + ytaxf - 12 (Aa^ + yi,axf 

3 

- 4 (Aa^ + t/^oa)^ + 3251^ ^ (^777| + ml). + tti^, + 2M| j - 2Kia^M^ 

2=1 

3 

+ 125(2 ^ ^ 4” '^‘hv- 4“ ^5 4“ ‘^M^'j + — 2Xi(i\^M2 

i=l 
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4 
— ( 
5^ 


3 r 


:gl ^ {^'nd + + m| + 2Mf j + 2/s;^ ^m| + + 2Ml^ 

i=l ^ 

+ 3 (a^, — 2\ia\^Ml^ + 2 (a^, — 2KiaKiMi^ 

+ ^ 2A^ (Qi + Q 2 ~ QV) + ^5 + 2Mf 

— 2Qs\\ (^QlTT^%d + Q2 ^‘h" + Qs^^s^ + 3^^ (Qd + Q^ ~ Q‘s) 




X + 2Mf j - SQs^^ (Qs^s + Qo^nj^^ + 

+ 2 (oa; ~ 2AiaAiM() (Qi + Q 2 ~ Ql) + 3 ~ 2K,iaf^^M[) (Q]j + — Q|) 

- ‘^^Qsg'iVt iQ2m]j^ + Qqml^ + Quml^) - 2AQsgiyl {Qim]j^ + Qgm^^ + Qdm\) 

3 

- l/r {Qi^k + <5L"ii3 + + ?>2Qsglgi ^ + Quml. 

+ Qdml. + + Qo^k) + ^‘^Qsghi 

3 

+ Y1 + Qimljd + Q2m 

2=1 
3 

+ ^ ^ (2QdTn‘^. + 2Qjjfn‘^ + 2QDinlj^^ + QQeTfi^^ + SQiTn^d + 3Q2Triffu + SQltti^ 


i=l 


Qw^w + Qh'^h, 


,2 


+ 1 3QjjTmjjr + SQH'm'H, 


2=1 


+ <5Q"^Qi + 8<5n"i«J 


+ 3(5s'iVf^^ + 3(5e + SQi + 

+ 2(5^ + 2(5^/ + SQi + 18Qg + SQI’ + 9(5^) + 2,Q^m^ + 2Q^,rri^/ 

3 

+ (3Qdmd, + ^Qk^k + ^Qk^k + '5e"ie, + “^Qlnikd + 2Qlm]ju + 2Qimi, 

2=1 

+ ^Qg'^Qi + Qs'^'si + ^Qu'^ui) “*" ‘^QsQk'^k' ‘^QsQk'^'n' 

3 

+ Qs (3Qdm^, + + 3Qkm% + Qlml. + 2Q\m]jd + 2g2rn|„ + 2Q2 


2=1 


+ ^Qq'^Q, + Qs'^Si + ^Ql'^Ui) 


(B12b) 


= 


i=l 


3 3 




+ mi .) + 4A, al, + ( mi. + + m^^ ) + 6/t, 


'i "-A 


‘Di 




'-i 


+ 8 EEA?^n + m\d + m^u + 2m|j + (AjOA^ + ^ja\kj 


i=l j=l 
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3 3 

+ 24 ^ ^ + ni]ju + 2m| + m|,, + + [Xia^. + Kjax,)^ 

i=i j=i 

3 3 

+ 18 ^ ^ [k-(^ 2m| + ml^ + m|, + + (/tia^. + Kja^y 

i=i j=i 

+ 12 A 2|/2 {2mjj^ + + m| + 771^3 + 771 ^ 3 ) + UX'^y^ ( 777 ^^ + 2ml^ + m\ 

+777^3 + 777^ + \X^yl {2mjj^ + mjj^ + ml + ml^ + mjj + 12 (Aa* + ytaxf 

3 

+ 12 {Xab + Vbaxf + 4 (Aa^ + yraxf - S2gl ^ (ml + m^^ + m^^ + 2M|j 

i=l 

3 

+®Ki ~ 27CjaKiM3] — 125f2 |^A^ 4“ ^77“ + + 2M|^ + — 2XiaXiM2 

i=l 

4 ^ 

— gfi'i ^ ^ |^3A^ (ml^d + 777^u + 777^' + 2iVf^^ + 3 c7^^ — QXiCiXiMi 

i=l 

+ 2k^ ^777| + 777|,^ + 777^^ + 2M'^ + 2a^, - AKitti^Mi 

+ 1^ — 4A^ (Qi + Q 2 + QI) + ^77“ 4“ 4" 2M(^ j 

i=l 

4- 2A^ (Qi + Q 2 + Qs) (Qs^^d + Qs^'h^ + Qs^'s 4" ^ 1 ^ ~ 4 (Qi + Q 2 4- Q|) 

^ (®Ai ~ 2AjaAiMQ + 2Qsa\. (Qi + Q 2 + Qs) ~ Stc^ (QI + Qh + 

X ^777| + 777|)^ + 777^, + 2M(^^ + StC^ (Qs + Qd + Qj)) (QS'fT^^S 4" Qs^\)^ + Qs^Jj^ + 

— 6 (QI + Qd + {al. — 27CjaK,M() + SQsal. {Qs + Qd + Qu) 

+ ^‘^Qsg'iVt {Q 2 + Qq + Qu) {^‘h^ + ^Qa + ^us) 4" ^‘^Qsg'iO^l {Q 2 + Qq + Qu) 

+ ^‘^Qsg'iVb {Qi 4" Qq 4- Qd) (^77d 4" + ^^ 3 ) 4- 12 ( 55 ( 71 ^ 0 ^ ((5i + Qq + Qd) 

+ ‘^Qsg'ivl {Qi + Ql + Qe) (^77d 4" + ^ 63 ) 4- AQsg'iol {Qi + Ql + (5e) 

— 96Qsglg'i {2 Qq + (5u + Qd + Qd + Qd) ~ 12(55(72fi'f 41^2 (9 Qq + ^Ql + 3Qi 

+ 3 Q 2 + Qjjt + Q 77 ') ~ gQsfi'iS'i^ 2Mi (2(5d + 2(5 d + 2(5d + 6Qe + 3Qi + 3 Q 2 + Qjfr 

+ Qh' + 3 Q 5 + Qq + 8 Q,) - SiSj] - 4Qsg[* [2Mf (QQ^ + 9Q| + 9Ql + 3Ql + 6Ql 

+ 6 Q 2 + 2Q^ + 2Qff, + QQl + 18Qq + 3Q| + 9Q^) + ( 6 Q 5 M(^ — S'l) Sq . (B12c) 

If the one-loop contributions to the effective potential from top and stop loops are also 
included, it is necessary to consider the expansions for yt, at, mg^ and 777 ^ 3 . The required 
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expressions for yt read 


d;’ = 


Vt 


+ %t +yl —~ ~ ^fi'i ~ ‘^9i {Q 2 + Q% + QV) 


13 


(B13a) 


=yt\^- 22i/i^ - ^yt - ^yhb - ylyl - + ^y'i + ^yl + yl + 2Sa + 3S^) 

+ 25'f [A^ (Q? — Q2 + Ql) + ‘^y't (2Qq + Qn) + l/b (Ql ~ Qq + Qd)] 

+ i%3i/t + QaM + al + ^yi^ + ^^3 + y^2 + 

+ 2fi'l‘^ [2 [Q 2 + Qq + Q^) + (Q 2 + Qq + Ql) ^q] + 8fl'3fi'2 + -^yial 


3913 

~450 


-9i 

136 


32 


fi'sfi'f {Qq + Qn) + fi'2fl'i + ^ 929 i {QI + Qq) 


^2/2 


I 16 

3Q^ + -QJ + -Ql + ( 3 Q 2 + Qq - 4QQ 


3'^^ ' 3 


(B13b) 


= yt\^^yt + yl/b + 13|/ti/fe + ylvl + Qa^ + Ibyl + byl + yl + 2Ea + 3E 
(A^ + ‘^yl + 12|/i) — 65f2 (A^ + l/fe + %t) — S'! 


16 2 
-y1?3 


52 2 


■ [A^ (Qi + 2Q2 + Q| + Qq + Q^) + yl (Qi + Q2 + 2Qq + Q^ + Q^) 


+12S? (Qi + Q^ + Qi)] + 1^91 - f 9l 

+ 2j;* (Qi + Qj + Ql) (Ql + Ql + Ql- Eg) + lesiji + ^dsl 


143 




208 


32 


ala'i {Ql + Qq + Ql) + -ralal + Qala'i {Ql + Qq + Q«) 


+ yE^^ 9? {Ql + Qq + Ql) |> 

and those for at read 
= «t 


(B13c) 


A^ + %t +yb ^al — ^al ~ ^^al — ‘^a'l (Q 2 + Qq + Q«) 

32 2 2 26 2 

+ yt I 2AaA + Vlytat + ‘2,yh(ib + ~^93^3 + ^92^‘i + 

+ ^9i {QI + Qq + Ql) , 

= atj - 22yt - 5yt - by^yl - ylyl - A^ (A^ + 3|/2 + Ayl + yl + 2Ea + 3E.) 
+ 2.gi [A^ {QI — QI + Ql) + 2yl (2Qq + Ql) + yl (Q^ — Qq + Q^)] 


13 


(B14a) 


..22 .22 2/62 2 2 \ 128 4 33 4 

+ 165'3l/t + %2l/t + 9 i ( gl/t + + yfi'2 + 


3913 


-9i 


46 



+ ‘^9i [2 {Q 2 + Qq + Qt) + {Ql + Qq + Ql) 

+ {Qq + Qu) + dldi + 65'2fi'f (Q 2 + Qq) 



SSy^at - 20ylab - lOytyb {ybat + ytab) 


~ 2ybyT {ybdr + yrdb) ~ 2\a\ (2A^ + 3y^ + 4?/^ + + 2 Sa + 

- 2A^ {3yt(it + 4:ybab + yrO-r + 2SaA + + 32g‘^yt [at — ytM^) 

+ 12fi'2l/t (ot - ytM 2 ) + ‘^gl [Qytat + 2ybab - {6y^ + 2yl) Mi] 

+ ^9i {Q\ ~ Q2 + QI) (®a — \M[) + 2|/t (2(5 q + (at — l/t^O 

^12 7H2fi 

+ (Q? -Ql + Ql) {ab - ybM[) ] - —gtM^ - QQglM^ - -^Q^Mi 

- Sg'iM] [2 {Q\ + Qq + Q'lOj + {QI + Qq + Q^) Sq] — 16(73(72 (1^3 + ^ 2 ) 

272 64 

—(^3 + (Qq + Qn) (^3 + M[) — 2glg‘l {M 2 + Mi) 

- I2glg'^ {Ql + Q^ (M2 + M^) - ^glg'^ (QQ^ + q| + leQ^) + M^ }, (B14b) 

32 26 

= Vt 2\a\ + 12ytat + 2ybab + —g^M^ + QglM 2 + —(73M1 
+ 4(7f M( (Q2 + Q^ + gy X A' + 67/2 + yl - ^gl - Sgj - ^gl 

- 2^f (Q2 + Ql + Ql) + at^5Ayf + y7/^ + ISy^yl + ylyl 

+ + ‘^Vb + 127/2) 

- 6(722 (A2 + 7/2 + 67/2) - 7/2 (|a2 + ^-yl + I772) 

~ 277(2 (A2 {QI + 2Q2 _|_ Q| _|_ _|_ _|_ y"^ {Q"^ _|_ Q2 _|_ _|_ Q2 _|_ 

+^‘2yl {Ql + Qq + Ql)] + ~^9t —^92 — 

+ 25'!^ (Q2 + Qq + Ql) (Ql + Qq + QI~ ^q) + ^^9^92 + ~^939i 
+ ■^5'35'f (Q2 + Qq + + ~^929i + 65'2fi'f (Q2 + Qq + l?n) 

26 1 r 

+ y^9i9i {QI + + Ql) r + l/( 1447/^0* + 2Aylab + I'iytyb {yt^b + o,tyb) 


+ 2ybyr {yb^T + abyQ + ISXyt {Xat + a^yt) + 8Xyb {Xab + axyb) + 2A7/^ (Aa^ + axyQ 


47 



+ SA^Oa + 4A ^ Ai (Aoa, + axK) + 6A ^ (Aa^^ + axKi) - GAgjyt {at - ytM^) 

i=l i=l 

32 

- -^ahb {ab - VbMs) - 36glyt {at - i/jMs) - Qgjyb {ab - i/bMs) - Qg^X {ax - XM 2 ) 

- ^alvt {at - ytMi) - ^glvb {ab - ybMi) - A {ax - AMi) 

~ {QI + Qq + Qu) (at — ytM[) — Ag'^yb [QI + Qq + Qa) {ab — ybM[) 

~ A (Qi + Q 2 + Q|) {a\ — XM[) + A8g2M2 + -^^g^Mi 

+ 8g[^M[J:Q (Ql + QI + QI)]. (B14c) 


The one- and two-loop /3 functions and the resulting O(t^) coefficient for toq^ are 


= ‘2y‘i {mjj^ + +ml^)+ ‘2aj + ‘^yl + 2^6 

<33 

- ^^glMl - QglMi - ^glMl - 8Qlg'^M[^ + + 2QQg'^i:\, (B15a) 

+ ml, +mX)- 2 ^yt {ml^ + ml^ + 771%) 

V3 

- ‘^ylyl (2m^^ + 771^3 + ml^ + mX + mlj - A^yja] - AQylal 

- 2 {yi,ar + yrabf - 2 X^y‘l (2m^^ + + m| + ml^ + 771^3) 

- 2A^7/6 {m]j^ + 2m]j^ + ml + 777^3 + 777 ^ 3 ) - 2 (Aof + ytaxf - 2 {Xab + yhaxf 

+ ^yWt {mX^ + 3"ig3 + ) + ^9i {at - ‘^ytatMi) 

+ ‘^ghb + ^Q3 + 4 M 2 ) + ^^2 (^^2 _ 2 y^abMi) + ^ ^2^ _ 3^^ 

+ H (""A “ ""S,) + 

i=l 

+ ((^2 - Qq + QX + ^Q3 + mX + 2Mf) 

+ 4c/f (Q 2 - Qq + Q«) («? - 27/tatM() - 24QQC/f t/^ (g2777^^ + Qqm^ + QumX) 

+ 45 'fl/fe (Ql — Qq + Qd) {mX^ + 777 q 3 -h 777^3 -h 2Mf) 

+ 4c/f (Ql - Qq + Qd) (flb - 27/bafeM() - 2AQQg'^yl {QimX^ + QgmX + QdmX) 

- ^Qgg'iyl {QimXa + QLmX + Qe^y 


4QQ^f YX (Qsml + QomX^ + Qd^|J + ‘^X^i 

i=\ 

3 

lOMi + ^ (2777^^ + 777^^ + 777^^ + mX^ + 777|J 


H: 


d + Q2mX^ + 


16 


-h 82glgl [M^ + M 2 M 3 -F M|) 
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ahl 2 {Ml + Ml M 3 + Ml) + 3 ^ - 2nil. - ml^ + m 


32^ 2 /2 

2 Qq939i 


f + M(M 3 + M|) + 3 ^ {2QQmQ. + + Qd^% 


if 3 

+ + 3^2 29M2 + tvM + ^ + ml. + m\i + 

-I _ ^=1 

1 r 3 

+ Talal 2 (Ml + Ml M 2 + Ml) + 3m|7 - 3m|, + 3 X] (^g, - - m%A 


+ ^‘^Qqala'i 2Qq (M(^ + M[M2 + M^) + Qjjrmlp + 

3 1 1 r 

+ ^ ( 3 gQm|, + Qiml^ + Qim]jd + Q2mjjn^ + —gt 289 Mi + 12 m |7 - 6m|, 

i=l 

3 -| 

+ ^ — 2ml). + 6m^, + 42mg, — 6m^d + 12m^u — 6m^, + 2mQ, — 24m^, j 

i=l 

+ Y^9i9'i ‘^QQ{^Qd + 9(52j — 9(5d + 9(5e ~ 9(5i + 9(^2 + ^Qw ~ ^Qh' ~ ^Ql 

+ lOQg — 18Qu) (Mi + Ml Ml + M^) + SQlpml^ — “^Qhimhi 
3 

+ 3 + Ql^% - Qlm^d + <52"^^“ 

i=l 

- Qlml^ - Qq^q, - + ^^QqQHrrriy + ^QqQmmjj, 

3 

+ Qq ^ {l2Qdml. + 12Q^m‘^^ + 2AQeml. + SQim^^a + l5Q2mjjn + SQLml. 

i=l 

+ 15QqmQ. + 12m^J + SQgfi'i'^ 3(5QMf (9(5^ + 9(5^ + 9Ql) + 3Ql + 6Ql 

+ QQ 2 + ‘^Qj[t + ‘^Qh' + SQi + 20Qq + 3(51’ + 9(5u) + 2Qljjrri^^ + 2Ql[,ml[i 
3 

+ (3(5d^dj + 3Qi)^i), + Qe^e; + 2Q\mljd + 2Q\mlju + 2Qlml. 

i=l 

+ 6(5Q"i|, + + 3(5n^nJ + ‘^QqQ^mM + ^QqQjj^mjj, 


Qq Y1 i^Qd^di + + ^Qo^Di + Ql^li + ‘2Qlm]jd + 2g2m^„ 


+ 2gimi^ + QQ^ml^ + Qlml. + 3g2i„y , 


(B15b) 


+ 771^3 + 771 ^ 3 ) + 2Ayl {m]j^ + 777^3 + 777 ^ 3 ) + A3ylal + A3ylal 
<33 


49 



+ ^yhl + ^%) + 4 {ytah + yhatf 

+ 2|/bl/r + ^q3 + + ml^ + ml^) + 2 (l/rOfe + ybarf 

+ ‘^>^'^yt + 2m|^ + m| + m|3 + ml^) + 2 (Aaj + ytaxf 

+ 2 A 2 |/fe ( 2 m^^ + mjj^ + m| + 771^3 + 777^3) + 2 (Aa;, + T/^aA)^ 

- + 2^3) - y£/3 («? - ‘^ytatM^) 

- Y^hb + "^d3 + 2 ^ 3 ) - y t/3 (ofe - ‘^ybabMs) 

- Qglyt + "^ q 3 + + 2 ^ 2 ) - 6^2 («? - ‘^ytatM 2 ) 

- Qghl {m]j^ + 777 ^3 + 777^3 + 2 Ml) - Qgl {al - 2 ybabM 2 ) 

- Y^ht + "^Q3 + + 2^1) - («? - 27/iatMi) 

- + "^Q3 + + 2 ^n - {(4 - ‘^ybabMi) 

- ^gh\ [Ql + + gy (rnk + rnl^ + ^^3 + 2 Mf) - 4^f (g^ + Q J + g^ 

X (a^ — 2ytatM['j + 2g'^y1 {Q 2 + Qq + Qu) {QQq'bn-jj^ + QQq'bn-Q^ + QQqw^^ + S'^) 

+ 12 gQ 5 (fa^ (g 2 + Qq + Qu) — ^g'lyl (Qi + Qq + QX {^‘na + ^Qs + ^X + 2 Mf) 

~ 45ff {Qi + Qq + Qd) (®6 ~ ‘^ybdbM'^ + 2g'^yl {Qi + gg + g^) 

X (6gQ777^^ + QQQmQ^ + QQqm^^ + S'J + 12QQg'^al {Qi + Qq + Qd) 

+ ^Qqg'iyl {Qi + Ql + Qe) + '^X) ^Qqgi^ (Qi + Ql + Qe) 

3 

+ ‘^Qqg'i |^ 2 A^ (gi + g 2 + Qs) ^ 77 “ ^s) 2®Ai (Ql + Q 2 + Q 5 ) 

^=1 

+ (Qs + Qd + Qu) {^4 + ^Di + (Qs + Qs* + Qd) 

- QQQqglg'i Ml {2 Qq + Qu + Qd + Qd + Qu) ~ 

96 

- ^2QQglg'^Ml (9Qq + 3Ql + 3Qi + 3 Q 2 + Qu7 + Qu) + ^S'l (^1 ~ ^ 1 ) 

~ [l2gQiVf^(2Qd + 2Qjy + 2gsi + 6Qe + 3Qi + 3 Q 2 + Qju “I" Qt^' + 3Ql 
+ Qq + 8Qu) + AM^ (3Q^ + 3Q|j — 3Q^ + 3Qg — 3Q^ + 3 Q 2 + Q^ ~ QX' ~ 3Q^ 

+ 3Q^ - QQl) - (6QqSi + s;) Sj] - 4gQ[2Mf (gQ^ + 9Q^ + 9Q1, + 3Ql + QQl 
+ 6 Q 2 + 2Q^ + 2Q^, + QQl + 18Qq + 3Q| + 9Q(^) + ( 6 QqM(^ — S'J Sq] . (B15c) 
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Finally, the relevant expressions for the soft mass read 


J + Aa^ - ^glM^ - - SQlg'^M'^ 


32 


32 


15^ 


— + ‘^Qug'i^'i, 


(B16a) 


/ 3 ll^ = -^ 2 yt {mjj^ + ml^ + mlj - Ay^yl + 2771^3 + ml^ + m%) 

- - 4 hjtab + Vhatf - 4:X‘^yt (2m|^ + m]j^ + m| + 171 ^^ + ml^) 

- 4 {Xat + ytaxf + l^glyl + tuq^ + ml^ + 2Ml) + 12^2 («? “ 2ytatM2) 

- ^ght (-5^k - ^Q3 + + 2^1) - {at - 2ytatMi) 

- \ 9 lyl (3m^^ - 771^3 - 2777^ - \glyl {mj,^ + m% - 2ml) 


8 2 [ 2 

-of y kX 


■ 3i [ 


n, m 


‘'Di 


777 


D, 


5^ 


+ y2/ 2 


H? 


777 


m 


2=1 


+ ^9iyt {QI + Qq ~ Ql) {'o^Hu + ^Qs + ^U3 + 2Mf) 

+ %f (<32 + <3 q - <3n) («t - 2ytatM[) - 2AQugiyt {Q 2 ml^ + Qq^Iq^ + (3n^^3) 

- 24Qug'M {Qi^Ha + <3Q"iQ3 + Qdml) - 8Qug'hl {Qi^]ii + <3L"ii3 + <3err7y 

3 

- 12(3n^f X] (<3s"i| + <3 d^d, + <3d^|,) 


2 = 1 
3 


3 

64 


^3 


2 2 

£<3^1 


- 8(3„^f ^ {Qi^Hd + <32^i„ + Qsml) 

' ' lOMi + + ml) 

2=1 

3 

8 (M^ + Ml M 3 + M|) - 3 ^ (ml - m^. + m^^ + m^^ - 2m 

2=1 

16g„ (Mf + M^Ms + Ml) 

I + Qoml^ + Qd^I. + ^Qqml. + Quml 

3 

I (2 2,2 2 

+ 

2=1 

3 

Qh'^^H' ^ ^ + 3QqT?2, 


2/2 

2Q22fi'3fi'l 


+ 12 E 


2=1 


12 

y 


2 2 
929i 


mjj, - m^p - 


+ ^2Quglg'l 
4 


2=1 


75 


£<i 


1261M^ + 21777^, + 3777 


2 

77' 
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3 1 

+ ^ + 4m^ — 12mg, + 21m^d + + 21m|^ + 3mQ, + 64m^, j 

2=1 

+ ~ 8Qu{9Qd + QQtj — 9(5d + QQe — QQi + 9 Q 2 + 3(5/^ — 3(5_h'' — 9Ql 

+ 9(5q ~ 22(5«) (^1 + MiM^ + + 12Q‘^,m^fji — 12Q^m^ 

3 

+ 12 ^ ( - - g^m|^ + g^m^^ - g^m^, + Qlmla - QWhu + gimi^ 

2=1 

3 

- Qq^Qi + Ql'^l) - ^^QwQnm^ + S3QH'Qum% + QuJ2i~ ^^Qd^di 

2 = 1 

- + 30g^m|,^ - QQeml. + 33gim^d - 15Q2mjju + 33QLm\. 

- 2igQjnQ. + 84g„m^J + ^QuQi 3QuM'^(9Q‘^ + 9Q^ + 9Qjj + 3gg + GQf 

+ QQ 2 4“ ‘^Qjp 2g|f/ + sg^ + isg^ + 3g^ + iig^) + 2Q^/tti^i + 2g^?Ti|p- 

3 

+ {^Qd^di + + ^Qn^Di + Ql^X + + 2 Qlmjju + 2gim^, 

2=1 

+ ^Qq'^^Qi + Qs^'si + ^Qu^ui) + ‘2QuQ^h'^‘h' + ‘^QuQ‘^rn^ 

3 

+ (Sg^m^, + 3Q‘^m^^ + SQ^ml. + Qlml. + 2Qlm]jd + 2Qlm]ju + 2gimi^ 

2 = 1 

+ ^Qq^q, + Qs^Si + , (B16b) 


= 48|/,^ + m ^3 + ml^) + Aylyl {m]j^ + m]j^ + 2ml^ + ml^ + m%) + 961 / 4 ^ 0 ,^ 

+ 4 (i/tOfe + + AX^yl {2m]j^ + + m| + 771^3 + + 4 (Aoj + T/tOA)^ 

- ^gWt {^‘hu + "^Qa + + 2^1) - y 1/3 («? - ‘^VtatMa) 

- ^‘^glyt + w?Q^ + m ^3 + 2Ml) - 12gl {a^ - 2ytatM2) 

- ^ght + "^Qa + + 2^1) - y l/i («? - 2i/iatMi) 

~ Sfi'fl/t (Q 2 + Qq + + ^Qa + '^U 3 + 2Mf) 

- 85ff (g2 + g^ + Qu) {o-t ~ ‘^yta-tM[) 

+ ^g'lyt (.Q 2 + Qq + Qu) {^Qu^nu^ + ^QuItiq^ + ^QuTri^^ + S'^) 

+ 12gnfi'fai (g 2 + gg + g«) + 32Qugiyl {Qi + Qq + g^) + ^^^ 3 ) 

+ ^‘^Qug'iOjb (Qi + Qq + Q^) 
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+ {Ql + Ql + Qe) ^Qug'lQ^ {Ql + Ql + Qe) 

3 r 

+ ‘^QuQi {Qs + Qd + Qij) j + Sa^i (Qs + Qo + Qd) 

i=l ^ 

+ 2A^ (Ql + Q 2 + Qs) "A “1“ ^ 5 ^ “I" Q"^ Qs) 


— M| (2Qq + Q^j + Qd + Qd + Qu) 

384 

— 12Qug\gi M 2 (9Qq + 3Ql + 3Qi + 3 Q 2 + Qh' + Qjd) — 2 ^^^ ^Mf) 

+ '^dlg'i (3Qd + 3Q^ ~ 3Q|, + 3Qg — 3Q^ + 3 Q 2 + Q|v — Qh' ~ ^Q1 

+ 3Qg — QQu) + ^QuM‘^{ 2Qd + 2Qj2 + 2Q£) + 6Qe + 3Qi + 3 Q 2 + Qw Qw 
+ 3Ql + Qq + 8QQ + (3Q„Si - 2S;) Sg] - AQ^g'^ [2Mf (QQ^ + ggs^ + gg^ 
+ 3Qg + QQ\ + 6 Q 2 4" 2 Q^ + 2Q^/ + 6 Q^ + 18Qq + 3Q|' + 9Q^) 


+ (6Q„Mf-S;) Eg], 


(B16c) 
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